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Thermodynamik des ebenen Ising-Onsager-Dipolgitters +) 
A.M. Dycune und J. B. Rumer 


Die Thermodynamik des ebenen Dipolgitters ist von groBem Interesse; denn dies 
ist das einzige Modell eines Systems, das einen Phaseniibergang zweiter Art er- 
leidet, der eine exakte Lésung zulaB8t. Die Lésung des Problems wurde in den 
Arbeiten [7] und [2] gegeben. Die in diesen Arbeiten entwickelte Methodik laBt 
aber prinzipiell keine Verallgemeinerung auf den Fall des dreidimensionalen 
Gitters zu. 

Eine Ubersicht tiber die genannten Arbeiten gab einer der Verfasser dieses Aufsatzes 
(J. R.) im Jahre 1954 [3]. 

Vor kurzem erschien eine Arbeit von Hurst und GREEN [4], die groBes Interesse 
beansprucht und wesentlich neue Gedanken fiir das Herangehen an die Lésung 
dieser Aufgabe enthalt. Die neue Methode ergibt eine einfachere Lésung als [7] 
und [2], und — was die Hauptsache ist — es ist die Méglichkeit nicht ausgeschlos- 
sen, sie auf das dreidimensionale Gitter zu verallgemeinern. Daher hielten wir es 
fiir zweckmaBig, die neue Methode in Form einer Ubersicht darzustellen. 


§ 1. Die Zustandssumme 


Gegeben sei ein rechtwinkliges Gitter, das aus m Spalten und » Zeilen besteht und 
N=nm Gitterpunkte enthailt. Wir numerieren die Gitterpunkte in lexiko- 
graphischer Ordnung, d. h., wir bezeichnen jeden Gitterpunkt mit der Nummer 


k=(k,—1l)m+h3° lakoim, 124,<-n, (1) 


“wo k, die Nummer der Zeile ist und k, die Nummer der Spalte, in denen der be- 
treffende Gitterpunkt liegt. 
Um Randeffekte zu beseitigen, verwenden wir periodische Randwerte; d. h., wir 
schlieBen das Gitter in beiden Richtungen ab, indem wir die (m + 1)-te Spalte 
mit der 1. Spalte und die (n + 1)-te Zeile mit der ersten Zeile identifizieren. 
In jeden der Gitterpunkte setzen wir einen Dipol, dessen Achse senkrecht auf der 
Gitterebene steht. Jeder der Dipole kann eine der beiden méglichen, einander ent- 
gegengesetzten Orientierungen haben. Offensichtlich ist die Gesamtzahl der még- 
lichen Konfigurationen der Dipole in dem Gitter gleich 2”. 
Zur Beschreibung der verschiedenen Konfigurationen gehen wir folgendermafen 
vor: Wir schreiben jedem Dipol eine diskrete Veranderliche o zu, die nur zwei 
Werte annehmen kann. Es ist o = + 1, wenn der Dipol nach rechts orientiert ist, 
und o — — 1, wenn der Dipol nach links orientiert ist. 


1) A.M. Dycune i J. B. Rumer, Termodinamika ploskoj dipol’noj resetki Izinga-Onsagera. 
Erscheint gleichzeitig in den Uspechi fiz. Nauk. 
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Wenn wir nur die Wechselwirkung zwischen den am nachsten liegenden Nachbarn 
beriicksichtigen und die zwischen den weiter voneinander entfernten Dipolen ver- 
nachlassigen, kann die Konfigurationsenergie H(o,..., oy) in folgender Form 
dargestellt werden: 


k=N 
H (040% careiis Oy) ee (J On On+4 + Jo Op Onem): (2) 
Darin ist -- J, die Wechselwirkungsenergie zwischen zwei benachbarten Dipolen, 
die in ein und derselben Zeile liegen bei gleicher bzw. entgegengesetzter Orientie- 
rung der Dipole. Analog bezeichnen wir mit -+ J, die Wechselwirkungsenergie 


zwischen zwei benachbarten Dipolen, die in ein und derselben Spalte liegen. 
Fiir die Zustandssumme erhalten wir den Ausdruck*): 


Z(T) = > exp [— Ho... cy)/T] = 


k=N 
=>) exp {- 5 (9, 0% On+1 + Oy O% oxen = 


N 
= 2 Lee (— 0, 04,0441) EXP (— 9 0% 0% +m)- (3) 


Darin bedeutet SS die Summation iiber alle 2% moéglichen Konfigurationen der 


(9) 
Dipole. AuBerdem wurden die Bezeichnungen 
0, = i (4) 


eingefiihrt. Beriicksichtigen wir, daB oj = 1 ist, so erhalten wir die Formel 


1+ o,0,tha« 


exp (x 0; 0,) = cha: 0,0, sho = (i th? ay (5) 


Mit den Bezeichnungen 
w= tl My th 0, 


k6énnen wir Z(7') in folgender Form schreiben 


N 
(DY Uy ae ead ae a i — © On OK+1) (L —Y OL on+m). (6) 


§ 2. Die Methode der Graphen 
Das in (6) vorkommende Produkt 
k=N 
G (x, Y, 6... Gy) ede — % 0% 0444) (L — Y 0% On+in) (7) 


ist ein Polynom der Veranderlichen 2, y, 0,, ..., oy. Da G in (6) hinter dem 
Summenzeichen der Summe iiber alle Konfigurationen steht, liefern nicht alle 


1) Statt k7’ wird stets nur 7’ geschrieben (Anm. d. dtsch. Red.). 
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Glieder des Polynoms einen von Null verschiedenen Beitrag zur Zustandssumme. 
Wir werden zeigen, daB nur diejenigen Glieder einen von Null verschiedenen 
Beitrag geben, die alle Veranderlichen 0,,..., Oy in der 0., 2. und 4. Potenz ent- 
halten. Es geniigt offensichtlich, diese Behauptung fiir irgendeine ausgewiihlte 
Veranderliche zu beweisen, etwa fiir o;. 

Wir sondern in G alle Faktoren aus, die o; enthalten. Offensichtlich sind dies vier 
Faktoren, und wir kénnen (6) in der Form 


Z (2; y) = [(1 es a) (i mre PNP oe G® (aa Y, Oro o Chis OER Gad o o) x< 


Oy +++ Of_y Gj4 4° ON 


x 2, 6: 6; 0-1) (1 — #6; Oj41) CL — y 07 6;—-m) (Lim Of Opens (9) 


schreiben, wo wir mit G® die Gesamtheit aller Faktoren bezeichnen, die die Ver- 
anderliche o; nicht enthalten. Fiihren wir in (8) die innere Summation nach 0; aus, 


k +7 k k+7 
k+m-2 kem-1 
: ; k+m+7 ktm k+m47 
eeegete ee : eo) shel eee ys 22 A+2m-7 + 
k+2m-1 k+2m 
e ° e e e es Ji sniae PBs» 
Bild 1a : Bild 1b Bild 1¢ 


so iberzeugen wir uns davon, daB nur die Glieder des Polynoms G einen von Null 
verschiedenen Beitrag zu Z(x, y) liefern, die o; in der 0., 2. und 4. Potenz enthalten. 
Diese Uberlegung gilt offenbar fiir alle Veranderlichen o;. 

Da oj = 1 ist, liefert jedes Glied des Polynoms, das alle Veranderlichen oy, . ++) ON 
in geraden Potenzen enthalt, nach dem Summieren tier alle Konfigurationen 
einen Beitrag, welcher der Gesamtzahl der Konfigurationen 2% proportional ist. 
Wir kommen zur Methode der Graphen, indem wir bemerken, dai man jedem 
~ Glied des Polynoms G(x, y, 01, ..., oy) eindeutig eine Gesamtheit von Linien 
zuordnen kann, die gewisse Paare von benachbarten Gitterpunkten verbinden. 
So entsprechen z. B. den Graphen von Bild 1 (a, b, c) die Glieder des Polynoms: 


2 
la On Ok+4 Sktmt1> 
D2 opel B: 4 2 2 o 
1b Ok Ok+1 Ok+m—1 Ch+m Ch+m+1 Ck+2m—1 Ck+2m> 
Dap 2 2 2 2 2 2 2 
Ke Of Oh+4 Onea=3 Ok+m—1 Ck +m Ck+mt1 Ck+2m—2 %k+2m—-1 [8k +3m—2 9K+3m-1- 


Wie Formel (7) zeigt, entspricht jeder horizontalen Linie der Faktor x und jeder 
vertikalen der Faktor y. ee : 

Wir haben oben bewiesen, da nur diejenigen Glieder des Polynoms G einen von 
Null verschiedenen Beitrag zur Zustandssumme beisteuern, die die Veranderlichen 
61, -+-, Oy in der 0., 2. und 4. Potenz enthalten. a 

Fashisch bedeutet dies, daB aus der Gesamtheit der méglichen Graphen nur 
diejenigen einen von Null verschiedenen Beitrag zur Zustandssumme beisteuern, 
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bei denen in jedem Gitterpunkt 0, 2 oder 4 Linien zusammentreffen. Einfachste 
Beispiele solcher Graphen zeigen Bild 1b und 1e. Alle Graphen, die einen von Null 
verschiedenen Beitrag liefern, miissen also geschlossen sein ; dabei ist es in einigen 
Gitterpunkten médglich, daB sie sich einmal selbst schneiden (Gitterpunkt k + m 


von Bild 1b). 
Die Zustandssumme kann also in folgender Form dargestellt werden: 
F(a, s4) 2 [Mi ) ay tS) (9) 
S(x, Y) = 2, dap. yh (10) 


Darin ist g,p die Zahl der geschlossenen Graphen, die aus « horizontalen und 6 
vertikalen Linien bestehen; dabei wird bei der Zahlung jeder Graph mit vielen 
Verbindungen als einziger gezahlt (z. B. der Graph von Bild 1c). 


§ 3. Einfiihrung von Fermi-Operatoren 


Zur Berechnung der Summe S (a, y) nach Formel (10) fiihren wir fiir jeden Gitter- 


3 
punkt zwei Paar Fermi-Operatoren Ons a, und 6,, 6, ein, welche die Bedingung 


A + + 
Ay Ay + A Cy = din} by b; + 6; by = Oi (it) 


erfiillen. (Alle anderen Antikommutatoren sind gleich Null.) 
Wir wollen nun den Ausdruck 


* ee +, 
S@7) sp IT (1 + 4x y ays, Op+m Oy b, + 42y ay Oy + Qe Opt m + 


gh + 
4 Qa dy Dyim + 2Yy by Oe+4 He ee ae eo Ube by) (12) 


berechnen und beweisen, daB S(x, y) = S*(x, y) ist. 

Wir bemerken noch, da hier und tiberall spater der Bequemlichkeit halber mit 
dem Symbol Sp C die GréBe Sp C/Sp H# bezeichnet wird, wo E die Einheitsmatrix 
der Dimension C ist. Bei dieser Schreibweise wirkt sich die Dimension der hinter 
dem Spurzeichen auftretenden Matrizen nicht auf die Resultate aus. 

Wir lésen das hinter dem Spurzeichen stehende Produkt auf und fassen die Glieder 
mit gleichen Potenzen‘(2a)* (2y)’ zusammen. Dann erhalten wir 


S*(x, y) = 3’ ay uy)? Sp ¥ ie 3) 


Darin stehen hinter dem Spurzeichen Summen verschiedener Produkte der Fermi- 
Operatoren Qj, a, und b;, bj. 

Um die Struktur dieser Produkte iibersetzen zu kénnen, ist wieder die Verwendung 
von Graphen bequem. 

Zu diesem Zweck fihren wir die in Bild 2 dargestellten acht elementaren Graphen 
ein, die jedem der acht Summanden zugeordnet werden, aus denen die Faktoren 
in dem Produkt (12) bestehen. 

Wir wollen Linien, die von einem Gitterpunkt aus nach oben und links verlaufen 
immer ausgezogen zeichnen und ihnen Erzeugungsoperatoren zuordnen, Die 
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Linien, die vom Gitterpunkt aus nach unten und rechts verlaufen, zeichnen wir 
gestrichelt und ordnen ihnen Vernichtungsoperatoren zu. Ferner setzen wir zu 


jedem Erzeugungsoperator oy den Faktor 2% und zu jedem Erzeucu t b 
den Faktor 2 y. J zeugungsoperator b;, 


—— --- 
| 


xe 


rey + | 
4 4XY 7 b..m ax b, 4XY ay by, Dfa7 Sag 
k 
k Qe ——e——— k 
| 
+ bee 
+ b i A ovspod 
2x04 Of. my 29% Oh 7 Exh 7 A “kek 


Bild 2 


Wir ordnen jedem Gitterpunkt je irgendeinen der acht elementaren Graphen von 
Bild 2 zu. 

Wir sprechen von einem ,, Doppelgraph‘‘, wenn in ihm zwei beliebige benachbarte 
Gitterpunkte entweder iiberhaupt nicht miteinander verbunden oder wenn sie 
durch eine ausgezogene und eine punktierte Linie miteinander verbunden sind. 


+ 
oak | z (kK) (k-m) 
—_ t-- i ee 
. | 
| C 


k1) @ | i 
| a) k b, | 6, 
| 
| a) b) (kK) 
pentane ——t SSS 
| 
{ 
| b) | 
Bild 3 Bild 4 


Wir heben hervor, da8 wir iiberhaupt keine Graphen finden kénnen, in denen 
einige der Doppellinien beide ausgezogen oder beide punktiert sind; denn aus 
Bild 2 ist zu ersehen, daB in den elementaren Graphen alle ausgezogenen Linien 
nach der entgegengesetzten Seite gerichtet sind wie die entsprechenden gestrichel- 
ten Linien. 

Wir werden nun beweisen, daB alle Doppelgraphen geschlossen sind. Wie Bild 2 
erkennen la8t, geht namlich von jedem Gitterpunkt eine gerade Zahl von einfachen 
Linien aus. Diese Linien werden nur durch Linien verdoppelt, die von benach- 
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barten Gitterpunkten ausgehen. In den ‘Doppelgraphen lauft also durch jeden 
Gitterpunkt eine gerade Zahl von Doppellinien, d.h., ein solcher Graph ist 
geschlossen. ) ‘ 
Offensichtlich entspricht jedem Doppelgraph ein gewisses System von. gepaarten 
Erzeugungs- und Vernichtungsoperatoren. Zum Beispiel kann die horizontale 
Doppellinie, die den (& — 1)-ten mit dem k-ten Gitterpunkt verbindet, konstruiert 
werden, wenn wir dem (k — 1)-ten Gitterpunkt einen der Graphen von Bild 3a 
zuordnen und den k-ten Gitterpunkt einen der Graphen von Bild 3b. 
Den uns interessierenden Teil des sich ergebenden Graphen zeigt Bild 4a. 
Analog kann eine vertikale Doppellinie konstruiert werden, wenn man dem 
(k — m)-ten Gitterpunkt einen der Graphen von Bild 5a zuordnet und dem k-ten 
Gitterpunkt einen der Graphen von 
Bild 5b. Der Graph, der sich dabei er- 
| | gibt, ist in Bild 4b dargestellt. 
=F iat ; Offensichtlich ist das konstruierte 
1 1 System der Doppelgraphen dem System 
der Graphen des vorhergehenden Para- 
graphen vollig aquivalent, und es geht 
in dieses ttber, wenn man jede Doppel- 
linie durch eine einfache ersetzt. 
sal | | AuBer den bis jetzt betrachteten Dop- 
—- —— pelgraphen wird es auch Graphen geben, 
| ; in denen einige von zwei benachbarten 
J b) Gitterpunkten durch eine einfache (aus- 
Bild 5 gezogene oder gestrichelte) Linie mit- 
einander verbunden sind. 
Augenscheinlich entsprechen solchen Graphen Produkte von Fermi-Operatoren, 
die nicht gepaarte Operatoren (Erzeugungs- oder Vernichtungsoperatoren) ent- 
halten. 
Setzt man solche Produkte aber in die Formel (13) hinter dem Spurzeichen ein, so 
ergeben sie zu der Summe den Beitrag Null, und man braucht sie daher nicht zu 
beachten und nur diejenigen Produkte in der Summe zu belassen, denen ge- 
schlossene Doppelgraphen entsprechen. 
Wir haben fiir S(x, y) und S*(x, y) die Ausdriicke 


S(x,y) = 2 Ia.p y?, (10) 


Wo gap die Zahl der geschlossenen Graphen ist, die aus « horizontalen und £ verti- 
kalen Linien bestehen, und 


S* (2, y) = > (22) (2y? Y Sp (dx... a,.:..8, canbe...) (14) 


or 


: : ; ea 

wo hinter dem Spurzeichen gepaarte Produkte von « Fermi-Operatoren a, und B 
+ 

Fermi-Operatoren b,, stehen, 


Wegen der bewiesenen Aquivalenz der einfachen und der Doppelgraphen wird die 
Zahl der Glieder in der inneren Summe genau g, B- 

Beriicksichtigen wir, da& bei der Bildung der Spur von Produkten aus Fermi- 
Operatoren gilt 


Sp Gy da = Sp A}, aa = 4, (14a) 


— 
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so erhalten wir nach Formel (14) 


pee 2 eye (ae LE. ae), (15) 


wo in der Klammer g, Glieder stehen, von denen jedes gleich +- 1 ist. 

Das Vorzeichen des betreffenden Gliedes ist dadurch bestimmt, ob man eine gerade 
oder ungerade Anzahl von Inversionen ausfithren muB, um von der urspriinglichen 
Anordnung der Fermi-Operatoren, die durch die Formel (12) bestimmt ist, zu der 
»gepaarten™ Anordnung zu kommen, bei der die Erzeugungsoperatoren neben den 
entsprechenden Vernichtungsoperatoren stehen. 

Im Anhang I wird bewiesen, da8 bei der Paarung der Operatoren das Vorzeichen 
des Produkts immer erhalten bleibt und da8 infolgedessen in der Formel (15) 
tiberall das Plusvorzeichen zu nehmen ist. Damit ist bewiesen, da8 


S(x,y) = S*(a, y). (15a) 


§ 4. Berechnung der Funktion S (a, y) 


Zur Berechnung der Funktion S(x, y) nach Formel (12) miissen wir die Spur eines 
Produkts von NV Faktoren berechnen, von denen jeder eine Summe von 8 Gliedern 
ist: Kins, eines Produkts von vier Fermi-Operatoren und sechs Produkten von zwei 
Fermi-Operatoren. 

Im Anhang II wird eine Formel abgeleitet, nach der man bequem die Spur eines 
Produkts von Faktoren berechnen kann, die Fermi-Operatoren linear enthalten, 
d. h. von Produkten der Art 


i=N 
IT Py (ain Ce + dix “| . (16) 
i=1 [® 


Wir k6énnten diese Formel benutzen, wenn es uns gelinge, jeden der Faktoren 
in (12) als Produkt von Faktoren darzustellen, die in bezug auf die in ihnen vor- 


kommenden Fermi-Operatoren linear sind. 
Um dieses Programm der ,,Linearisierung“ des Produkts in (12) durchzufihren, 


_berechnen wir das Produkt 
(AP + Qed) (AP + 2y8,) (AP + ares) (AP + Berm), (17) 


P : = 4 4 
wo Ay vier Operatoren sind, die mit den Operatoren az, a, by, dx antikommu- 
tieren. Wir erhalten 


+ + 2 
AP AP AP AD + day dy by ays1 Opem + A? AR Qp+1 Oetm + 
ne 
+ AP AP 2ybeem by + AWAY 2y br bevy + 


(4) 


cs + 
+ AS A® 2a Oy Dem + AP AP 22 Gy 1 ay + 4AP AP xy db. (18) 


. Oa : 
Wenn wir den vier Operatoren A\ die Bedingungen 


Sp (AP AP AP AM) = 1; Sp (4p Ay) =1; 4u¢=1,2,3,4 (19) 
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auferlegen, so konnen wir S(x,y) in der Form eines Produkts von 4N Faktoren 
darstellen 


k=N A & : 
S(«, y) =Sp [] (A® + 22&) (AP + 2ydy) (AP + ess) (AP + Prem). (20) 
k=1 


Die Bedingung (19) kann man erfiillen — wie leicht nachzupriifen ist —, 


: Dae 
wenn man sechs weitere Paare von Fermi-Operatoren Psd é ee (A = 1, 2, 3, 4, 5, 6) 
einfiihrt und setzt: 


AP = 2 (co + of? + cf”) 
AP = 72 (oP + of? + o¥”) 
Ap = 8 (6 + a + of) et) 
Ap = 2G +e + 2). 
Mit den Formeln (20) und (21) kann man die Spur des Produkts (13) als Spur eines 


> 
Produkts von 4.N-Faktoren darstellen, die in den Operatoren ay» ay, by, Ox, Cee 


cy, linear sind. Nach Formel (12) von Anhang II kann das Quadrat dieser Spur 
als Determinante einer 4N-reihigen antisymmetrischen Matrix D geschrieben 
werden, deren Elemente durch den Doppelindex (k, 4) numeriert werden, wobei 


k=—2,1,..., Nund/d = 1, 2,3, 4 ist. Die Elemente der Matrix D sind: 
(EAD A) = 1 (A6  A) (ep Ane in A a ed ee 
(b,3|D\k+1,1)=2 (k+1,1/D\/b3)=—e (22) 


(k,4|D| k+m,2)=y (k+m,2|D|k, 4) = —y. 


Die iibrigen Matrixelemente sind gleich Null. Wir sehen, daB die Matrix D als 
eine Ubermatrix dargestellt werden kann, deren Elemente vierreihige Matrizen 
sind: 


EAE ey word TATE 
XT “iyiraeee Ory a0 
eed Meee eee he Mais (23) 


Darin ist nach den Formeln (22) 


ah ee 00 x0 0000 

ro Daan ate Weg ['0.208.05.0)\a8 ePID OORN 
“Tea OF ke oe aL es | (24) 

1 0000 0000 

—1 -1 -1i 0 0000 0000 


Die Struktur der Matrix D zeigt, daB sie als Summe von direkten Produkten 


N -reihiger Matrizen EZ, E, und E,, jeweils mit den vierreihigen Matrizen U, X, Y 
geschrieben werden kann 


D=EXU+E,XX+E, x ¥Y—EP x X?— ET x YT, (26a) 
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wobei gilt: 
(lB) k)= dn; @|ALE+1)=44; ((|Enlk+m)= 6. (25) 
Nach Anhang II (12) ist 
S*(e. 1) t= Det DD) (26) 


Bei der Berechnung der Determinante der Matrix D kénnen wir sie einer Ahnlich- 
keitstransformation unterwerfen, die Formeln von Anhang III verwenden 


D’ =(S x E) D(S2 x #) =EXU+M x X+ BAY — EP XT EP TYT (27) 


und die Matrix S so wahlen, daB ZH, und E,, Diagonalmatrizen werden. 
Wenn wir bedenken, daB N = m-n eine groBe Zahl ist, so kénnen wir den 
Elementen der Matrix D die Periodizitatsbedingung 


WY se b A| DIN -- kW) = (bh, A | DE, 2) (28) 

auferlegen. Dann gilt 
Ey, = (H,)", B= (B)Y; (B44 = Et = Ef (29) 
(E,)¥—" = Bo” = BE. 


Hieraus folgt, daB die Eigenwerte der Matrizen HZ, und E,, Wurzeln N-ten Grades 
aus Kins sein werden: 


Qt 


(P 1E| p) = exp ("5") = a9: (p [Ey |p) = exp (= mp) = Br. (30) 


Die Matrix D’ schreibt sich in Form einer Stufenmatrix 


Uy 
Dice |Maaicy aye} (31) 
‘Us 
wo nach (27) und (30) gilt 
* 
U,=U + aX + f,¥ —a,X? — p,¥* (32) 
oder nach den Formeln (24) 
0 1° 1+ a, 1 
_ 0 1 1+ 6,4 : 
a a= 1 Boy (33) 
* —1l—ajx —-—2 0 1 
1 1 —Bpy 1 0 


Berechnen wir die Determinante, so erhalten wir 


2 


2 es 2 ad — 4) cos — mp 
eat) (li ya 2 y (1 — 27) cos =e p ae y a : 
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Wir fiihren die Bezeichnungen ein 


Pp =Pan+ Dy 1< py; Op, Sn 


2 27 p 7 ea 
C08 =, —p = COs - (». | AF cos W mp = oo a Py 


(34) 


und schreiben den abschlieBenden Ausdruck S?(a, y) in Form eines zweifachen 
Produkts : 


1=2 P2=™M Dy) 
(S(z, y))? = Th t + a2) (1 + y2) —20 (1 — y?) cos = p, — 
Pi=1p,=0 n 
— 2y (1 — x?) cos ss @ — alk (35) 


Durch Logarithmieren erhalten wir hieraus 


Ps=mM 


—1 22 
In a 4 a) (1 + y2) — 2x (1 — y*) cos ~ p, 
p.=0 


—2y (1 — 28) e0s = “lp, a alt 


In S(x,y) = e py 


(36) 
Bei groBen n, m geht diese Summe in ein Integral tiber 
In S(x,y) = saan [ [(t + 2%) (1 + y2) —22 (1 — 9?) cos oy — 
— 2y (1 — 2?) cos w,] da,, dw, (37) 


§ 5. Die freie Energie; der Phasentibergang 


Da wir einen Ausdruck fiir In S (x, y) haben, kénnen wir unter Verwendung von 

Formel (9) den Ausdruck fiir die freie Energie (Ff (a, y) in der Form schreiben 
F (x, y) 

TN aa N In Z (x, y) = 


1 
In2 — In [(t — 24) (1 — )] + 


2 an 


1 + 2) (1 + y?) —2a (1 — y?) cos a, — 


— 2y (1 — 2*) cos w,] dw, da, (38) 
Setzen wir x = th 6, und y = th 8, ein, so erhalten wir 
F(e,y) _ rp 


1 
— ay In 2 4 5 af fm [cos 26, ch 26, — sh 26, cos w, — 
0 0 


—sh 26, cos w] dw, dag 


(39) 
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genet Ausdruck wurde in der urspriinglichen Arbeit von ONSAGER [1] abge- 
e1tet. 


Wir wollen nun beweisen, da der Ausdruck, der in (38) hinter dem Logarithmus- 
zeichen auftritt, nicht negativ ist. 
Es ist 


[(1 — 2?) (1 — y*) — 4ay]? = 


= (1 = 2)? (1 — 93)? + 160% 42 — Say (1 — 2%) (1 — 42) > 0 (40) 
oder 


(1 + at — 22%) (1 + yt — 242) > Bay (1 — 2%) (1 — 2) — 160% y?. 
Addieren wir auf beiden Seiten der Ungleichung den positiven Ausdruck 
4a? (1+ y!) 4 4? (1 + at) = 4a? (1 — y2)? 4 4y? (1 — 22)? + 1622y2, 

so erhalten = 

(1+ at + 20%) (1 + yt + 29%) > (2 a(t —y®) + Dy (1 —29)) 
Hieraus folgt 

A (x, y) = (1 + 2*) (1 + y’) — 2a (1 — y*) — 2y (1 — 2?) 2 0. (41) 
Dabei wird das Gleichheitszeichen bei den Werten x = x, und y = y, erreicht, 
fiir die gilt 

(1 — x¢) (1 — y¢) —4ay. = 9. (42) 

Verwenden wir (41), so kénnen wir die freie Energie in folgender Form schreiben: 


F (x,y) 
TN 


= R (x,y) + 


| In [A (ay) + 2% (1 — y?) (1 — cos 4) + 
0.0 


+ 2y (1 — 22) (1 — cos ,)] da, dos, (43) 


i a! 
2 (22) 


a 


wo R (x, y) eine fiir x = 2, y = y, regulare Funktion ist. 

Wir sehen: Eine besondere Eigenschaft der freien Energie hangt damit zusammen, 
daB der Ausdruck unter dem Logarithmus Null wird fiir @, = @, = 0 und bei der 
Temperatur 7',, die durch die Formel 


i 
th =a ht =y, (44) 


bestimmt ist, wo x,, y, die Beziehung (42) erfiillen. 
Da der Ausdruck 


A (x, y) +2a (1 — y®) (1 — cos @,) + 2y (1 — x?) (1 — COS wg) 
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fiir die Werte w, = 0, 2 =0,4= %, y= Ye ein Minimum hat und Null wird, k6nnen 


wir ihn in eine Reihe entwickeln und erhalten: 


alee A} cou +2(so52 )ae aay + 3 (Ay) J +a, (1 — y2) wf + ye (1 — 2) w3- 
(45) 
Beachten wir, daB nach (44) 


2) 1 2) 2 
de=—(1— 28) 47: dy = —0 — Ry Gp 


ist, so erhalten wir fiir die freie Energie in der Umgebung des singularen Punktes 


2a 2a 
P (2) pate 
— AD = Rl, w) + x5 spay | | meen r= re 
ir (Lye) Oe ae Yo Le re) ws] dw, dor, (46) 


1 024 024 : ; 
Cr.) — 272 E Gale — ae) — 2d, df (args). — a) (1 = Y:) + 
a2 : 
(5G) ct — a2 (47) 


Integrieren wir (46), so erhalten wir fiir die freie Energie den Ausdruck 


F 
Na Ry (PO) ae Gy ey i ae (48) 


Fir die innere Energie U erhalten wir nach der Formel 


U 1 F 
eee (x7) = Ry (T,) + C,(P.) (1 —T,)In(T—7,)? (49) 


und bemerken, daf sie im Punkt 7’ = 7’, stetig bleibt. Die Warmekapazitat er- 
leidet aber einen logarithmischen Sprung 


6 d {U 
N dt 


al = Ry (Tyee 0.) a eas (50) 


Schlu8 


Zum. Schluf bemerken wir, daB die geschilderte Methodik nicht unmittelbar auf 
den Fall eines dreidimensionalen Gitters verallgemeinert werden kann. Wir 
numerieren die entsprechenden Formeln fiir das dreidimensionale Gitter mit den 


gleichen Zahlen wie im zweidimensionalen Fall und versehen sie mit einem 
Stern. 
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Wir geben jedem Gitterpunkt die Nummer 


k = (kz — 1) mn + (ky —1)m+ ky Le ren 
1<k<n (1*) 
N = mnp es keep. 


Wir kommen leicht zu den Formeln 
Z (x, y, 2) = 2% [(1 — a?) (1 — y?) (1 — 2)]-4? 8 (a, y, 2), (9*) 


Se, 9,.2) = Bi Gany mya. (10*) 
a,P,Y 


Darin ist gag, die Zahl der geschlossenen raéumlichen Graphen mit «, f,y Linien in 
drei aufeinander senkrechten Richtungen. . 


Weiter kann man analog zu (12) Fermi-Operatoren einfiihren, die drei aufeinander 


senkrechten Richtungen in dem Gitter a;, a,, by, bi, Om. Cf entsprechen. Die elemen- 
taren Graphen analog denen von Bild 2 kénnen auch 

leicht konstruiert werden. Man kann sich leicht da- is) 6 
von tiberzeugen, da es im raumlichen Falle 32 ele- 

mentare Graphen gibt. Fiir die GréBe S* kann eine 

Graphenmethode entwickelt werden, wobei sich leicht 7 

die Formel 


Mere yet ee? (a 1 tt) 1). (15*) 

a,B,y & 
ergibt; dabei stehen in der Klammer g,,, Glieder. 
Zum Unterschied von dem ebenen Fall kann man 3 4 
aber hier nicht beweisen, daB alle Summanden in Bild 6 
(15*) positiv sind. AuBerdem gibt es ein einfaches 
Beispiel, das zeigt, daB dies nicht der Fall ist. Hierzu dient der Graph von Bild 8 


+ + + se ee ea ee 
* 
OG, Oz be ln, 0,07 CO, 0g Cy 070, 0a Ca. (=) 


Fihren wir die Paarbildung so durch, wie es (*) angibt, so erhalten wir, daB der 
Graph von Bild 6 einen negativen Beitrag zu der Zustandssumme beisteuert. Diese 
Tatsache hangt nicht von der Reihenfolge der Operatoren in den Eckpunkten ab, 
da ein Eckpunkt jedes Typs in dem betrachteten Graphen zweimal vorkommt. 


Anhang I 


Wir wollen beweisen, da8 bei der angenommenen Ordnung der Operatoren in (12) 
alle Graphen einen nicht negativen Beitrag zur statistischen Summe beitragen. 

Um von einem beliebigen Glied der Summe in (14) zu den Graphen tibergehen zu 
konnen, mu man die Fermi-Operatoren so umstellen, das wir von der lexiko- 
graphischen Anordnung, die in (12) angenommen wurde, zu einer Ordnung 
kommen, wo jeder Erzeugungsoperator neben dem entsprechenden Vernichtungs- 
operator steht. Wir bemerken, daB die Anordnung, in der die entsprechenden 
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Erzeugungs- und Vernichtungsoperatoren hintereinander stehen, wegen der 
geltenden Beziehung F 

Sp ay, ay = Sp Ay Ay 
unwesentlich ist. ' a 
Es geniigt, den Beweis fiir zusammenhangende Graphen zu fihren, da ein nicht- 
zusammenhangender Graph gleich dem Produkt der zu ihm gehorenden zu- 
sammenhangenden Teile ist. Der Beweis wird durch Induktion gefihrt. 
Hierzu iiberzeugen wir uns davon, dab 
der einfachste Graph, den Bild 1b zeigt, 
das Pluszeichen hat. Diesem Graphen 
entsprechen die Operatoren 


4 
Opty Ox+m On+m+1 Ut 1Uet+m1 
e 


+ $e + 
Ditm Ue+mts OK+m+i- 


Fiihren wir die Umstellung in der ge- 
wohnlichen Weise aus, so kénnen wir 
uns davon iiberzeugen, daB der Graph 
das Vorzeichen plus hat. 

Wir betrachten zuerst einen beliebigen 
zusammenhangenden Graphen, der sich 
nicht selbst schneidet. Offensichtlich 
kann ein beliebiger derartiger Graph 
aus einem elementaren durch aufein- 
anderfolgende Anwendung einer der 
Operationen von Bild 7 konstruiert 


apes 
Cte 
C-- S9 
> 


SS a werden. 
Es genitigt nachzuweisen, da eine jede 
d) derartige Operation das Vorzeichen des 


Graphen nicht andert. Da der elemen- 
tare Graph in (15) einen positiven Bei- 
trag liefert, zeigen wir damit, daB alle 
Graphen, die sich nicht selbst schneiden, 


Bild 8 


einen positiven Beitrag liefern. Der Beweis wird fiir alle Operationen von Bild 7 
analog geftihrt. Daher betrachten wir hier irgendein Beispiel, etwa 7e. 
Wir schreiben die Operatoren heraus, die dem linken Graph von 7e entsprechen. 


+ + + 
0 Ot Opim B OKs m A+m+1 U+m+2U+mt+1 Y- (1) 
SSS a el 
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Darin bedeutet « die Gesamtheit der Fermi-Operatoren, die zu den Gitterpunkten 
gehéren, die dem k-ten vorangehen, # diejenigen, die zu den Gitterpunkten 
zwischen dem k-ten und dem (& + m)-ten gehéren (in Bild 7e ist die betreffende 
Stelle des Graphen durch eine Wellenlinie bezeichnet), und y die zu Gitterpunkten, 
die auf den (k + m + 1)-ten folgen. 

Nach Hinzufiigen des elementaren Rechtecks geht der linke Graph von Bild 7e 
in den rechten tiber. Wir schreiben die entsprechenden Operatoren heraus 


++ = 
X Ag+, Ay Opty Opin bk tmad Ak+m+2)- ; (2) 


a, : ng ; + 
Paaren wir in (2) die Operatoren y+, byim SOWIC @y+m+y> U¢+m+ 1, SO erhalten 
wir 


sign (1) = — sign « ay B Cee y. (2a) 


Nach demselben Verfahren iiberzeugen wir uns, daB 


sign (2) = — sign a dy B dg+m+1 Y- (2b) 
Somit erhalten wir 
sign (1) = sign (2), (2c) 


und die Graphen von Bild 7e haben das gleiche Vorzeichen. 

Auf dieselbe Weise wird der Beweis in allen anderen Fallen gefiihrt. 

Um den Beweis fiir die Graphen zu fiihren, die sich selbst schneiden, mtissen wir 
in die Zahl der auszufiihrenden Operationen auch Operationen vom Typ des 
Bildes 8 einbeziehen. Dabei wird der Beweis genau analog dem vorhergehenden 


Fall gefihrt. 


Anhang Il 


Es mégen 2t-Formen gegeben sein, die in den Fermi-Operatoren c¢,, c, linear 
sind, 
k=N 
L; 


I 


(ain Cy + Bin cr): (3) 


Wir leiten eine Formel zur Berechnung der Spur eines Produkts von derartigen 
linearen Formen ab 


sie — ts OURS Oa aeae PE (4) 
Aus der Definition (3) folgt: 
[Lj Lie = Dy Ly + LiL; = > (aye dive [Cx ev] + bie Caw (eu ew') = 
eh 


= >} (aj bin + G7 djx)- (5) 
E 


Sp LL; =4 > (ain bin + Gin bjx) = 4 (Li Lil. 
3 
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Wir stellen die Folge der 27 — 1-Identitaten zusammen: 


Sp (Ly L,L3..-L_7) + Sp (LL, Lg. -- Lex) = [L, L.] Sp (Lg Ly - - - Lex) 
—Sp (L,L,L,... Ly) — Sp (LpL3 Ly... Lo2) = —[L,L,] Sp (LL, « ae L,;) 


Si aug Get leriSiva) @lie vo 0.66.08) a) wes B Jolseng: ie] ie: si (6: @, oe fol (al w Koyo) eriegion-e veloc) ©: eltie7. 1 (¢hk@ ste ih%s SSS Retin aoa eae 


Be. e ple ahéue 0141 01.0) 88 6 8) 0 Sle le © © loi 16 1816) © celal ie. (0,6) Oe) euielie’ isp elieme tenon en eu ohne. SuSricike go Ci: A colo Salas ck coe eae 


Beriicksichtigen wir, daB sich die Spur des Produkts bei zyklischer Vertauschung 
der Faktoren nicht andert, so erhalten wir, indem wir (6) addieren, 


O'S (Ly Lae) = > (=I) [Lal Sin, aoe ee eens 


So, = De (—1)' Sp (Z, L;) Sp (LZ, 23... Le ee ee ae) (7) 


1=2 


eine Rekursionsformel zur Berechnung der uns interessierenden Spur. Wenden 
wir die Formel (7) nacheinander an, so erhalten wir 


Sor = 2" (—1)? Sp (Ly, Le.) - Sp (Le. Lu.) «Sp (Ler-1 Lar). (8) 


Darin bedeutet das Zeichen ’, daB sich die Summation itiber alle Permutationen 


ky ke 2 PR Kory Kor 9 
2 ee oe eeoe Cue. 
erstreckt, die die Bedingungen 
ky <i kgs+ hy <3 Kane ge vice ee eee (10) 


erfiillen, und P ist der Permutationsgrad von (9). 


Ein Ausdruck von der Form (8) wird als Pfaffscher Ausdruck [5] der dreieckigen 
Anordnung der Zahlen 


Sp (LZ, L,) Sp (Ly VEE Oasys (Ly Lg7) 


Sor == | iw Se BR ee eee (11) 


bezeichnet. Als Minor des Pfaffschen Ausdrucks, der dem Element Sp (LZ, L;); 
j > entspricht, wird der Pfaffsche Ausdruck bezeichnet, der sich aus den 
Pfaffschen Ausdruck (7) dadurch ergibt, daB8 wir die i-te und die j-te Zeile sowie 
die i-te und die j-te Spalte der Formel streichen. (7) ist die Entwicklung des 
Pfaffschen Ausdrucks (11) nach den Elementen der ersten Zeile. 

Das Fundamentaltheorem der Theorie des Pfaffschen Ausdrucks, das wir hier 
ohne Beweis anfiihren, stellt einen Zusammenhang zwischen dem Quadrat des 
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Pfaffschen Ausdrucks und der Determinante her, die der antisymmetrischen 
Matrix 
( 0 Sp (Z,L, Sp (L,L5) ... Sp (L, Lor) 
—Sp(LZ,L,) 0 Sp (ZL) ... Sp (LZ, Lor) 
(Spx)? = |— Sp (Z,L,) —Sp(Z,L,) 0 ... Sp (L_Ler) (12) 


SPS, 6) Sle: Olli Ol (60.0 6M nie! ee) 6) 8) 0 ie Bf pie) 0) 8! o” 0) 0. Se ens By 3.6) bifecete, 6 


entspricht. 


Dadurch ist die Berechnung der Spur nach Formel (4) nach den Formeln (5) und 
(12) auf die Berechnung einer Determinante zuriickgefiihrt. 


Anhang III 


Das direkte Produkt einer n-reihigen Matrix A mit der m-reihigen Matrix B ist 
die n+ m-reihige Matrix: ; 


LIE [cate Ce aed v3) 
ALB 2.38 AyjywB 
Offensichtlich gilt A x B=: Bx A. 
Aus der Definition folgen die Formeln fiir die Produkte: 
(ATSB) (DseDY= AC SCBE 
CAS BG sD). (FIG) = ACK Oe BDG:, (14) 


Das direkte Produkt kann man einer Ahnlichkeitstransformation unterwerfen 
(Sie 8.) (A 5) (Site a) =. 8,45," 2C SBS, . _ (15) 
In dem Spezialfall S, = H# erhalten wir die Formel 
(BERS) C4 SOB) (ES 8) =A XS ,BSy, (16) 


die im Text angefthrt ist. 
Ubersetzt von F. BarrEts 


Deutsche Red. durch D. FRanK 
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§ 1 Einfithrung 


In den letzten Jahren hat sich der Anwendungsbereich der Magnetohydrodynamik 
stark erweitert. Man findet sie heute in der Theorie der mit fliissigen Metallen 
gektihlten Kernreaktoren [7], der Theorie der magnetohydrodynamischen Me8- 
gerate [2], der Sonnenphysik [3], der Geophysik und Astrophysik [4], der 
Theorie der Gasentladungen und des Plasmas [5—7], der Theorie der gesteuerten 
thermonuklearen Reaktionen [8] und der Magnetoaerodynamik [9, 10]. 

Die Gleichungen der Magnetohydrodynamik sind nichtlinear. Das fithrt zum Auf- 
treten einer Reihe spezifischer Effekte, insbesondere zur Bildung von StoBwellen. 
Diese StoBwellen spielen eine wesentliche Rolle bei der Kompression und Auf- 
heizung eines Plasmas [1/—13]. Sie stellen einen der Mechanismen zur Bildung 
interstellarer Magnetfelder dar [74]. Ohne StoBwellen ware es unméglich, die Uber- 
schallumstrémung eines Kérpers im Magnetfeld zu behandeln [15, 16]. Magneto- 
hydrodynamische StoBwellen sind auch wesentlich in der Theorie der Impulsent- 
ladungen im Plasma [77]. 

Im vorliegenden Artikel wird, bei Horrman und TELLER [78] beginnend, ein Uber- 
blick uber eine Reihe theoretischer Arbeiten gegeben, die den magnetohydro- 
dynamischen StoBwellen gewidmet sind. 


§ 2 Einfache Wellen 


Eine wichtige Klasse nichtlinearer Lésungen der magnetohydrodynamischen 
Grundgleichungen bilden die einfachen Wellen. Als einfache Wellen bezeichnet 
man solche Lésungen, in denen alle magnetohydrodynamischen GréBen, so die 


1) Uspechi fiz. Nauk 72, 33 (1960). 
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Dichte 9, der Druck p, die Entropie s, die Materiegeschwindigkeit v und das 
Magnetfeld H, als Funktionen einer einzigen (z. B. o, das selbst nur von der Ko- 
ordinate x und der Zeit ¢ abhangt) auftreten. Wir beschranken uns auf den Fall 
ebener eindimensionaler einfacher Wellen. Die besondere Rolle, die einfache 
Wellen in der Magnetohydrodynamik spielen, besteht darin, da® beim Fehlen 
eines Sprunges nur sie an das Gebiet einer Gleichstrémung angrenzen kénnen 
[19, 20]*). Kine spezielle Art von einfachen Wellen sind die Automodellwellen, 
d.h. Wellen, bei denen die magnetohydrodynamischen Gré8en nur vom Verhalt- 
nis x/¢ abhangen. Automodellwellen entstehen immer dann, wenn die Anfangs- 
bedingungen von Parametern frei sind, die die Dimension einer Linge haben. 

In der Magnetohydrodynamik gibt es drei Arten von einfachen Wellen [19, 20]: 
1. Alfvénsche Wellen 


cH 
y4 0 : 


He z= Vm + HZ = const, 9 =const, s=const, H, = const. 


C= y= yw + vz = const, (2.1) 


Es gilt « = + 1 fiir Wellen, die sich in positiver z-Richtung ausbreiten, e = — 1 
fiir Wellen, die sich in entgegengesetzter Richtung ausbreiten. 
2. Magnetoakustische Wellen 


dye sl db; H, , U, 
=e—, Se = (2.2) 
do e” de 4x? (U2, — Uz) 
cH, U2H, dp 


Py As w) d =c*, s=const, H, = const. 
e QU — Ve @ 


- Dabei ist U, die Alfvén-Geschwindigkeit, die bestimmt wird durch 
— Ht . 
/420 
v, und H; sind die Transversalkomponenten der Geschwindigkeit und des Magnet- 
-feldes, c ist die Schallgeschwindigkeit und 


yi 


Us, = V4 (U2 ++ (U2 + oP — 40203); 


die Zeichen Plus bzw. Minus bei der Gré8e U.. entsprechen einer schnellen bzw. 
~ Jangsamen magnetoakustischen Welle. 
3. Entropie-Wellen 


do = (00/0s),ds, p =const, v=const, H = const. (2.3) 


1) Die Behauptung von Muersasast [21], er habe eine Lésung gefunden, die foe os 
Welle ist, aber an das Gebiet einer Gleichstromung angrenzt, ist fehlerhaft. Damit die ewe- 
gung des Mediums sich so verhilt, wie sie Muarpasast angibt, muff man an den ay einen 
veranderlichen (und nicht einen konstanten) duBeren Druck anlegen, der sich nach einem 
speziell vorgegebenen Gesetz andert. Das wurde z. B. in einer Arbeit von A. G. KULIKOVSKIS 


[22] festgestellt. 
38* 
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Bei der einfachen Alfvénschen Welle sind die GréBen p, 0, Vz, UV, H nicht verander- 
lich, und die Vektoren » und H drehen sich um die z-Achse um jeweils gleiche 
Winkel. Die Phasengeschwindigkeit bei der Ausbreitung Alfvénscher Wellen ist 
gleich v, -+ ¢U,. Da sich die letzte GroBe nicht andert, breitet sich die einfache 
Alfvénsche Welle ohne Veraénderung ihrer Form aus. Daraus geht im besonderen 
hervor, daB die Alfvénschen Wellen keine Automodellwellen sein kénnen. 

Die Phasengeschwindigkeit bei der Ausbreitung von magnetoakustischen Wellen 
ist gleich v, + ¢U. Wenn, wie gewéhnlich, die adiabatische Kompressibilitat 


mit wachsendem Druck fallt 
(oh 1 
= i OF 
‘eal ( Q 


so folgt aus den Differentialgleichungen (2.2) 
d 
— (v, + eU,) > 0. (2.4) 
de 


Die Ungleichung (2.4) bedeutet, daB bei einer Bewegung der magnetoakustischen 
Welle ihr Profil verzerrt wird. Die Steilheit der Welle wird in Dilatationsbereichen 
(do/dt < 0) geringer, und in Kompressionsbereichen (de/dt > 0) wird sie gréBer. 
‘Das fithrt schlieBlich zur Bildung von Spriingen (StoBwellen) in den Kompres- 
sionsbereichen!) [26—29]. Die Lésungen der Gleichungen fiir einfache magneto- 
akustische Wellen haben eine sehr komplizierte Gestalt?). Wir wollen uns nur mit 
einigen qualitativen Schlu8folgerungen beschaftigen, die unmittelbar aus den 
Differentialgleichungen ersichtlich sind. Vor allem stellen wir fest, daB der Druck 
sich im gleichen Sinne andert wie die Dichte; das transversale Magnetfeld H, 
andert sich in einer schnellen magnetoakustischen Welle im gleichen Sinne und in 
einer langsamen im entgegengesetzten Sinne. Einfache magnetoakustische Wellen — 
sind eben; wenn zu Beginn v, = 0 und H, = 0 ist, so werden diese Verhaltnisse 
beibehalten. Wenn sich eine Automodellwelle nach rechts ausbreitet (ec = + 1), 
und wenn die GréBen H, und H, gleiches Vorzeichen haben, so wird die Longi- 
tudinalgeschwindigkeit v, kleiner. Die Transversalgeschwindigkeit v, wachst in 
einer schnellen Automodellwelle an; in einer langsamen nimmt sie ab. 


§ 3 StoBwellen. Zemplenscher Satz 


Wie wir gesehen haben, werden Spriinge im Kompressionsgebiet einfacher Wellen 
infolge Verzerrung des Profils gebildet. An der Sprungflache gelten dann die 
Gesetze der Erhaltung der Masse, des Impulses und der Energie. AuBerdem bleibt 
auch die Transversalkomponente des elektrischen und die Normalkomponente des 
magnetischen Feldes erhalten [78, 30, 31]. Wir beschranken uns auf eine Betrach- 
tung von Spriingen konstanter Amplitude. 


1) Fir den Fall, daB das Magnetfeld senkrecht zur Ausbreitungsrichtung der Welle steht, 
wurden analoge Resultate von S. A. Kapnan und K. P. SransuKxovié [23], SmGRE [24] und 
Tantutt [25] gefunden. 


*) Wie Frrepricus (s. [29]) zeigte, sind die Gleichungen fiir einfache Wellen durch reine 
Quadraturen zu lésen. 
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Kine Klassifizierung aller méglichen Typen von magnetohydrodynamischen 
Spriingen wurde von 8. I. SYrovarsxry durchgefiihrt [32, 33]. 
Wie bekannt ist [34, 35], gilt in der gewéhnlichen Hydrodynamik bei Erfiillung 


der Bedingungen 
en 1 0 ds 0 
ap}, \ o EAU ap ie (3.1) 


der Zemplensche Satz, auf Grund dessen sich in der StoBwelle der Druck und 
die Dichte vergréBern: py > p,, 02 > 0,, mit anderen Worten: StoBwellen sind 
immer Kompressionswellen.!) 

L. D. Lanpav und E. M. Lrrstc [30] haben nachgewiesen, daB auch in der 
Magnetohydrodynamik bei Erfiillung der Bedingungen (3.1) StoBwellen mit kleiner 
Amplitude bei willkiirlicher Richtung des Magnetfeldes Kompressionswellen sind. 
Eine Behauptung, die der Umkehrung des Zemplenschen Satzes entspricht und 
darin besteht, daB in einem idealen Gas in der Kompressions-StoBwelle die 
Entropie anwachst, wurde yon Horman und TELLER bewiesen [18]. 

Der Zemplensche Satz selbst wurde von Horman und TELLER nur fiir senkrechte 
StoBwellen mit kleiner Amplitude bewiesen; er ist aber bei willkiirlicher Intensitat 
der Welle und bei willkirlicher Richtung des Magnetfeldes [36 —38] giiltig. 


§ 4 Bedingungen fiir die Entwicklung der Spriinge 


Die Angabe der Grenzbedingungen auf dem Sprung reicht nicht aus, um eindeutig 
die Sprunglésung zu bestimmen. Dieser Schwierigkeit begegnet man auch in der 
gewohnlichen Hydrodynamik. So sind z. B. beim Herausziehen eines Kolbens aus 
einem Rohr formal Wellen méglich, die zwei Losungen entsprechen: 1. der Auto- 
modellverdiinnungswelle, 2. der VerditinnungsstoBwelle [39]. Die zweite Lésung 
wird in der Hydrodynamik abgelehnt, weil sie dem Zemplenschen Satz wider- 
spricht. In der Magnetohydrodynamik sind VerdiinnungsstoBwellen, wie schon 
weiter oben bemerkt wurde, ebenfalls unmoglich. 
Es gibt jedoch in der Magnetohydrodynamik eine groBe Anzahl von Kompressions- 
stoBwellen, und fiir das Problem, die Bewegung des Mediums bei vorgegebenen 
~ Anfangs- und Grenzbedingungen zu beschreiben, existieren unzahlige Losungen. 
Wenn man z. B. einen ideal leitenden Kolben in ein ruhendes magnetohydrody- 
namisches Medium hineinbewegt, in dem das Magnetfeld H, die Richtung der 
Kolbennormale hat, so sind in diesem Falle, wenn man die folgenden Ungleichungen 
aufstellt 

2 2 De Ler? 

EIEN Se ee AAO ea Tee ec, 
4U 0, /4U?2, — 3c 


(Fas 4s ol V4 709; Cy: Schallgeschwindigkeit ; der Index Null kennzeichnet das un- 
gestérte Medium; w: Geschwindigkeit des Kolbens; y = 5/3) zwei Typen von 
Lésungen méglich: 

1. die gleiche KompressionsstoBwelle wie beim Fehlen des Magnetfeldes, 


1) ZeMPLEN, C. R. Acad. Sci. (Paris) 141, 712 (1905); 142, 142 (1906). (Dt. Ubers.) 
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2. zwei magnetohydrodynamische StoBwellen, die sich mit gleicher Geschwindig- 
keit bewegen. Diese Geschwindigkeit hat die GroBe 


2 4u? 
Ue " rE / ge tee 
wobei das transversale Magnetfeld H,, zwischen diesen Wellen bestimmt wird 
durch die Formel 


/ Bxog( D217 ie ee 
Bh ae | ae 


Beide StoBwellen sind Kompressionswellen. Die Dichte des Mediums zwischen 
diesen Wellen ist gleich 


Die Entropie wachst in beiden StoBwellen an. 

So sehen wir, da8 die Entropiebedingung, mit deren Hilfe es gelang, in der ge- 
wohnlichen Hydrodynamik ,,iiberflissige’ Spriinge auszuschlieBen, sich in der 
Magnetohydrodynamik als viel zu schwach erweist. 

In Wirklichkeit konnen nicht alle StoBwellen, bei denen die Grenzbedingungen 
erfiillt sind und die Entropie wachst, zustande kommen. 

AuB8erdem ist namlich nétig, daB die Lésung stetig von den Anfangs- und Grenz- 
bedingungen abhangt, d.h., daB eine unendlich kleine Schwankung der magneto- 
hydrodynamischen Gro8en eine unendlich kleine Veranderung der Lésung her- 
vorruft. Solche Spriinge werden wir nach I. M. Gu’ranp [33a] evolutionare 
(evolucionnye) nennen. Bei nichtevolutionaren Spriingen ruft eine unendlich 
kleine Schwankung eine endliche Veranderung der Lésung hervor und auBerdem 
eine Aufspaltung des Ausgangssprunges in einige Spriinge von endlicher GréBe[40]. 
Nichtevolutionare Spriinge sind deshalb instabil in bezug auf die Aufspaltung und 
existieren nicht. Die Evolutionsbedingungen fiir die Magnetohydrodynamik 
wurden in [40] aufgestellt. Um festzustellen, ob eine magnetohydrodynamische 
StoBwelle evolutionar ist, mu8 man zu den Gleichgewichtswerten der magneto- 
hydrodynamischen GréBen p, 0, v,, Vy, v,, H,, H, die unendlich kleinen Schwan- 
kungen 0p, 09, dv;z, dvy, dv,, Hy, 6H, hinzufiigen, deren Abhangigkeit von den 
Koordinaten und der Zeit die Form 


exp 1 (ka — at) 


hat. Diese Schwankungen kann man als Superposition von Wellen unendlich 
kleiner Amplitude darstellen (Magnetoakustische, Alfvénsche und Entropie- 
Wellen). AuSerdem mu man noch die Wellen beriicksichtigen, die von der 
Sprungflache auslaufen. Nach der Linearisierung der Grenzbedingungen erhalt 
man sieben lineare homogene algebraische Gleichungen in bezug auf die Ampli- 
tuden der verschiedenen Wellen auf beiden Seiten der Sprungflache. Dabei muB 
man beriicksichtigen, daB die Geschwindigkeit der StoBwelle D auch den unendlich 
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kleinen Zusatz 6D erhalt. Nach Eliminierung der Gré8e 5D bleiben sechs unab- 
hangige Gleichungen (Grenzbedingungen). Wenn die Zahl der nach beiden Seiten 
von der Sprungflache auseinander laufenden Wellen auch gleich sechs ist, so er- 
halten wir ein System aus sechs linearen homogenen algebraischen Gleichungen _ 
mit sechs Unbekannten (den Amplituden der Wellen unendlich kleiner Intensitat). 
Voraussetzung fiir die Existenz einer nichttrivialen Lésung (in der die Amplituden 
der Wellen von Null verschieden sind) ist nun, daB die Determinante des Systems 
gleich Null ist. Die so erhaltene Gleichung stellt die Abhangigkeit zwischen @ und k 
dar. Wenn einem beliebigen & ein reelles @ entspricht, so ist der Sprung stabil im 
tiblichen Sinne; wenn aber méglichen Werten von k ein w mit positivem Imaginar- 
teil entspricht, so ist der Sprung instabil, die Anfangsschwankung wird exponen- 
tie mit der Zeit anwachsen. In beiden Fallen ist der Sprung aber evolutionar. 
Wenn auch die Schwankungen exponentiell anwachsen, so werden sie doch im 
Verlaufe einer kleinen Zeitspanne noch klein sein. 

Damit der Sprung nichtevolutionar wird, ist es nétig, daB die Stérungen, die im 
Anfangsmoment, t = 0 klein waren, zu einer beliebigen Zeit ¢ > 0 nicht kleiner 
werden. Hierbei ist die friiher durchgefiihrte Linearisierung der Gleichungen nicht 
zulassig, was dazu fiihrt, daB die Zah] der Gleichungen zur Bestimmung der Wellen- 
amplituden mit unendlich kleiner Intensitat nicht gleich der Zahl der Unbekann- 
ten, d. h. nicht gleich der Zahl der Wellen ist, die von der Sprungflache aus nach 
beiden Seiten auseinanderlaufen. Weil die Zahl der unabhaingigen Grenzbedin- 
gungen in der Magnetohydrodynamik gleich sechs ist, so besteht die Evolutions- 
bedingung darin, daB im allgemeinen die Zahl der auseinanderlaufenden Wellen 
auch gleich sechs sein mu8. In der Magnetohydrodynamik gibt es vierzehn ver- 
schiedene Phasengeschwindigkeiten fiir die Ausbreitung von Wellen mit unendlich 
kleiner Amplitude: 


I elcinag Jom ; 
aes ie ie oe se One Neve U2 — Onre Dia Mes 
Uy, — Uy-, ye, Vox + Ure, Vor — Vor, Vax + Uo,; Voz — Ur,, 
Veg + Us, Vor — Uns Vaz 


(der Index 1 bezieht sich auf das Gebiet vor der StoBwelle (x < 0), der Index 2 
auf das Gebiet hinter der StoBwelle (2 > 0); das Koordinatensystem wurde so 
gewahlt, daB in ihm der Sprung ruht und sich in der Ebene x = 0 befindet; die 
~ Richtung der a-Achse wurde so gelegt, da die Projektion der Geschwindigkeit des 
Mediums auf die x-Achse positiv wird). Den im Gebiet vor der StoBwelle ausein- 
~anderlaufenden Wellen entspricht eine negative Phasengeschwindigkeit, den im 
Gebiet hinter der StoBwelle laufenden eine positive Phasengeschwindigkeit. 

Von den vierzehn aufgezihlten Phasengeschwindigkeiten entsprechen vier Ge- 
schwindigkeiten einlaufenden Wellen: 


Ve + EE Gacate hess OF Che, Vig 
und vier gehéren zu auslaufenden Wellen: 
Vor She Oss Vox + ie Vor Sipe Vou 


(alle diese Phasengeschwindigkeiten sind offensichtlich positiv). Die tibrigen Wellen 
werden einlaufende oder auslaufende sein in Abhangigkeit von den Beziehungen 
zwischen den GréBen v, und U,, U,, U_. Auf Bild 1 ist die Gesamtzahl der aus- 
laufenden Wellen bei verschiedener GréBe von v,, und v,, dargestellt. Den evolu- 
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tiondren Wellen entsprechen die Gebiete der 7 ,-V2,— Flache, in denen die Zahl der 

auslaufenden Wellen gleich sechs ist. - 
Wie wir sagten, wird die StoBwelle evolutionaér, wenn das Problem der kleinen 
Stérungen eine eindeutige Lésung hat. Dazu ist notwendig, da8 die Zahl der Glei- 
chungen (Zahl der Grenzbedingungen minus eins) der Zahl der Unbekannten 
entspricht (Zahl der Amplituden aus- 

Vax laufender Wellen)*). 

Die Gleichheit der Anzahl] der Gleichun- 
poe eG Iee 206-8 | 24527 gen mit der Anzahl der Unbekannten 


Ure reicht jedoch fiir Existenz und KHin- 
SEED Mas Ohi eit Oeg aaneminge deutigkeit der Lésung nicht aus. Es 
Up iy kann vorkommen, da die Gleichungen, 
246 =81| 245 20a) 1+ J2Gu ated die zur Bestimmung von Amplituden 

U,_ auslaufender Wellen dienen, und die 
GHOST METS =O | IGGE OM Tape Grenzbedingungen in isolierte Gruppen 

y,, zerfallen. Hierbei miissen die Evolu- 


U;- Us Uy + tionsbedingungen (die Zahl der aus- 

Bild 1. Zahl der Wellen, die von der Sprungfliche laufenden Wellen sei gleich der Zahl 

der Alfvénschen Wellen, der zweitesimmand 4 Unabhangigen Grenzbedingungen) 

die Zahl der magnetoakustischen und En- nicht nur fiir die Gesamtzahl der Ver- 

tropiewellen, anderlichen, sondern auch ftir jede iso- 

. lierte Gruppe gesondert erfillt werden. 

Ein solcher Zerfall der Gleichungen und der Grenzbedingungen in zwei isolierte Grup- 

pen trifft in der Magnetohydrodynamik fiir Wellen zu, die sich senkrecht zur Sprung- 

flache ausbreiten*). Tatsachlich sind in Alfvénschen Wellen mit kleiner Amplitude 

die Werte dv, und 6H, von Null verschieden (wenn die xy-Ebene so liegt, daB H,=0 

ist); wahrend in magnetoakustischen und Entropiewellen die Werte de, 6p, 

OVz, Vy, 0H, von Null verschieden sind. In diese beiden Gruppen zerfallen auch 
die Grenzbedingungen, die linearisiert sind in bezug auf kleine Stoérungen: 


1. Alfvénsche Storungen 
He Of, 
yp 0) as * “= ’ 
ie 2 rp 0 
{v,6H, — H,6v,} = 0; (4.1) 


2. Magnetoakustische- und Entropiestérungen 


(g(6v~ — 6D) + v,69} = 0, 
{on + 2ov, (dv, — 6D) + v760 + am =i); 
TU 
H,6H 
OND. es ee 
{ers Vy ee 0, (4.2) 


{H,dvy — H,(dv, — 6D) — v,0H == 0; 
il 
0Vz [Vz (dv, — dD) + v,dvy + dw] + ves (v,H, —v,H,) on,) = (0 
( 6D: Storung der StoBwellengeschwindigkeit, w:Enthalpie). 
') In dieser Formulierung wurden die Evolutionsbedingungen yon Lax [19], K. I. Banenxo 


und I. M. Gex’ranp [41] angegeben. 
*) Dieser Umstand wurde zuerst von S. I. SyRovaTsKIs [42] bemerkt. 
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Die Grenzbedingungen (4.1) enthalten 5D nicht. Deshalb sind sie alle unabhangig, 
und die Zahl der auslauferiden Alfvénschen Wellen mu8 gleich zwei sein. Die 
Grenzbedingungen (4.2) enthalten die Stérung der StoBwellengeschwindigkeit 5D; 
nach ihrer Eliminierung verbleiben vier unabhangige Grenzbedingungen. Folglich 
mu die Zahl der auslaufenden magnetoakustischen und Entropie-Wellen gleich 
vier sein. Aus Bild 1 geht hervor, daB 


zwei Evolutionitatsgebiete fiir StoB- | 


wellen existieren!) (Bild 2): ex 
1. Schnelle StoBwellen (auf Bild 2 mit 
dem Plus-Zeichen gekennzeichnet) : 


(Oe << Via» Ose << Vox <= LR (4.3) 


2. Langsame StoBwellen (auf Bild 2 
mit dem Minus-Zeichen gekennzeich- 
net): 


: , T ; 

. - P UV; Ua U1 
Bis jetzt haben wit Evolutionsbedin- Bild 2. Evolutionsgebiete fiir Sto8wellen (schraffiert). 
gungen in bezug auf die Stérungen be- Das Zeichen + kennzeichnet die schnelle, das 

: mA Zeichen — die langsame StoBwelle. 
trachtet, die nur von vund t abhangen. Der evolutionére Teil der Stofadiabate wird 
Kine Berechnung der Stérungen von durch die ausgezogene Linie dargestellt, der nicht- 
alloemeine r - 2 evolutionare Teil ist gestrichelt wiedergegeben. 
ee Eiee A t (die auch yor Y und Mit dem Buchstaben4 wird der Alfvénsche Sprung 
abhangen), fiihrt zu den gleichen Evo- gekennzeichnet. 


lutionsbedingungen (4.3) und (4.4) [44]. 

Man mu8 betonen, daB ein wesentlicher Unterschied zwischen Nichtevolutionitat 
und Instabilitat besteht. Instabile Zusténde kénnen bei der Bewegung des magneto- 
hydrodynamischen Mediums unter der Wirkung innerer Ursachen entstehen. 
Sie existieren einige Zeit, bis die Fluktuation eine kritische Hohe erreicht hat, 
wonach der instabile Zustand zerfallt. Nichtevolutionare Spriinge kénnen nicht aus 
sich selbst entstehen. Sie konnen sich nur unter der Wirkung von auBeren Fak- 
toren (z. B. beim Zusammensto8 von Gasmassen) bilden und nur einen Augen- ° 
blick als Spriinge in den Anfangsbedingungen existieren, wonach sie schnell 
in einige StoB- bzw. Automodellwellen aufspalten. Kine solche Spaltung einer 
nichtevolutionaren magnetohydrodynamischen StoBwelle wurde in der Arbeit [45] 
_ untersucht (s. § 8). 

Mit der Methode, die in diesem Paragraphen dargelegt ist, kann man nach- 
weisen, daB auch alle anderen magnetohydrodynamischen Spriinge (Kontakt-, 
Tangential-, Alfvén-Spriinge) evolutionar sind. 


§ 5 Konsequenzen der Evolutionsbedingungen 


Die Bedingungen (4.3), (4.4) haben wichtige Konsequenzen. Wenn z. B. ZWwel 
StoBwellen gleicher Art (beide schnell oder beide langsam) aufeinander folgen, 
so holt die nachlaufende die vorlaufende ein. Zum Beweis dieser Behauptung 


1) Die Arbeiten [43] und [40] enthalten eine unrichtige Behauptung tiber die Existenz eines 
dritten Evolutionitatsgebietes. 

Oye <_< Uy4, Us-< ez < Uae: 
Dabei wurde aber der oben besprochene Zerfall der Grenzbedingungen in zwei isolierte Grup- 
pen nicht beachtet. (In den Arbeiten wird an Stelle von Evolutionitat der weniger zutreffende 
Terminus ,,Stabilitat in bezug auf Spaltung‘‘ benutzt.) 
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betrachten wir z. B. zwei langsame Wellen. Die Geschwindigkeit der vorlaufenden 
Welle ist in bezug auf das Medium, das zwischen den Wellen eingeschlossen ist, 
gleich v,,, und die Geschwindigkeit der nachlaufenden Welle ist gleich v,,. Wie 
aus der Ungleichung (4.4) hervorgeht, ist fiir Wellen des untersuchten Typs 
Vin > Uy, Ver < Uz-. Da sich die Geschwindigkeiten U,- und U,_ auf ein und 
dasselbe Gebiet beziehen, ist U; = U,— und folglich ry, > Vge- 

Bei Wellen verschiedenen Typs kann man leicht feststellen, da der Alfvénsche 
Sprung die langsame StoBwelle einholt und daB die schnelle StoBwelle alle Arten 


von Spriingen einholt. Analog kann man feststellen, daf eine StoBwelle einen — 


schwachen Sprung einholt, wenn er zum selben Typ gehért wie die StoBwelle, 
oder zu einem langsameren. Ein schwacher Sprung holt eine StoBwelle gleichen 
Typs und Sto8wellen langsameren Typs ein. 

Da ein schwacher Sprung eine Grenze zwischen einer einfachen Welle und dem 
Gebiet stetigen FlieBens darstellt, kommen wir zu der SchluBfolgerung, daB sich 
nach jeder Seite nicht mehr als drei StoBwellen oder einfache Wellen bewegen 
kénnen ohne einander einzuholen: vorlaufend eine schnelle Welle (StoBwelle oder 
einfache magnetoakustische Welle), hinter ihr eine Alfvénsche Welle (Sprung- 
oder einfache Welle) und schlieBlich eine langsame Welle (StoB- oder einfache 
Welle). 
Benutzt man die Evolutionsbedingungen und das Zemplensche Theorem, so 
kann man auch bestimmte SchluBfolgerungen iiber die Anderung des Magnet- 
feldes in schnellen und langsamen StoBwellen ziehen [38]. Wir benutzen die Be- 


ziehung 
H v2, — U? 
Lie = Des os my (5.1) 
Oia 02 Uie 


die eine Folge der Grenzbedingungen ist. Aus der Formel (5.1) folgt, daB sich bei 

schnellen StoBwellen das transversale Magnetfeld H, vergroBert und bei lang- 

samen verkleinert. In beiden Fallen andert das transversale Magnetfeld seine 
Richtung nicht?). Kleine Magnetfelder 

H 2 Ox DE 


ees 


werden beim Durchlaufen der StoBwelle verstarkt, waihrend groBe Magnetfelder 


Hy 5s. 01 Vie 
8a 2 
geschwacht werden. Das weist auf eine etwas ausgleichende Rolle der StoBwellen 
ae on Durchgang einer groBen Zahl unabhangiger StoBwellen tritt das 
a atische Gleichgewicht mit Gleichheit der magnetischen und der kinetischen 
nergie ein?): 
He Tgilas 
8a aes, 


(5.2) 


*) Die in den Arbeiten [46], [43], [47] enthaltenen Behauptungen iiber die Méglichkeit der 
Richtungsinderung des transversalen Magnetfeldes in der StoBwelle riihren von der Nicht- 
beachtung der Evolutionsbedingungen her. 

*) Die Relation (5.2) (Gesetz iiber die Gleichverteilung der Energie) ist charakteristisch fiir die 


Magnetohydrodynamik. Sie wurde, von anderen Uberlegungen ausgehend, in den Arbeiten 
[45—51] gefunden. 


—— 
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In der Literatur gibt es Hinweise dafiir, daB schwache Magnetfelder 
H? 
a <P (5.3) 


beim Durchlaufen einer StoBwelle verstarkt werden [14], [52] und [53], Seite 253. 
Genauer gesagt, existieren bei einer beliebigen Beziehung zwischen H, und p, 
zwei StoBwellen: eine schnelle, in der das Magnetfeld verstarkt wird, und eine 
langsame, in der das Magnetfeld geschwiicht wird. Jedoch darf bei Anwendung 
der Ungleichung (5.3) die langsame StoSwelle nur eine unendlich kleine Intensitiit 
haben, wie aus den Evolutionsbedingungen (4.4) folgt: 


0< —< ——. (5.4) 


Bei Verletzung der Ungleichung (5.4) hért die langsame StoBwelle auf, evolutionar 
zu sein und wird gespalten. Auf diese Weise kénnen in einem Medium, in dem der 
magnetische Druck bedeutend geringer ist als der hydrostatische (5.3), in Wirk- 
lichkeit nur schnelle StoBwellen existieren, bei deren Durchgang das Magnetfeld 
verstarkt wird. Aus dem Dargelegten folgt, da magnetohydrodynamische 
Wellen einer der Mechanismen zur Bildung interstellarer Magnetfelder sind [14]. 


§ 6 Die Stoadiabate 


In der Magnetohydrodynamik wird die StoBwelle charakterisiert durch die Werte 
aller magnetohydrodynamischen GréBen 01, p,, v;, H, vor der StoBwelle und 
durch den Wert einer von ihnen (z. B. 0,) binter der StoBwelle. Alle tibrigen GréBen 
hinter der StoBwelle, auch der Druck p,, sind Funktionen von a,. Die Kurve 
P= f (1/02) wird StoBadiabate genannt. Wenn man das Koordinatensystem so aus- 
wahlt, da8 in ihm der Sprung ruht und in der Flache 2 = 0 liegt, die Projektion 
des Magnetfeldes auf die z-Achse gleich Null ist und der Geschwindigkeitsvektor 
parallel zum Magnetfeld ist, so wird die Gleichung der StoBadiabate 


Als Parameter enthalt sie die GréBen p,, 0,, Hy, Ay: 
1 

wi (— Drona ales Hy). (6.1) 
02 


An Stelle der fiinf dimensionsbehafteten Parameter 05, p,, 0;, Hy, Hy, kann man 
zur Charakterisierung der StoBwelle drei dimensionslose Parameter benutzen [14], 
FEB AUC 0 yl Gre] Or wober, U die Alfvénsche Geschwindigkeit und Vig die 
Geschwindigkeit der StoBwelle in bezug auf die ruhende Fliissigkeit ist. 

Fiir das Auffinden der StoBadiabate ist es ndétig, in den Grenzbedingungen alle 
magnetohydrodynamischen GroBen zu eliminieren, die nicht in Formel (6.1) ein- 
gehen. Die Gleichung der StoBadiabate in der Magnetohydrodynamik hat folgen- 
des Aussehen [30, 31]: 


1 ¢ 
(e€. — &) 4 ptr a | ( ) (Ho, — Hy,)? = 9. (6.2) 
02 Q1 
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In dieser Formel wird die innere Energie e, durch p, und g2 mit Hilfe der Zustands- 

gleichung ausgedriickt; H,, setzt man entsprechend Formel (5.1) ein. 

Die Gleichung (6.2) ist eine algebraische Gleichung dritten Grades in bezug auf 

Px bzw. Q9; folglich gibt es drei Zweige der StoBadiabate [14, 46, 36, 54), jedoch 
nur zwei von ihnen sind evolutionar [55]. Im 


?, a) Spezialfall der senkrechten StoBwelle (H, = 0) 
| gibt es nur einen evolutionéren Zweig der 
Co ’ StoBadiabate. 

As we - C Die Anderung der magnetohydrodynamischen 
4 Ce eee GréBen in der StoBwelle wurde zum erstenmal 

oie von He.rer [14] untersucht. Die Darstellungen, 
l 
be ES 


3 79 7: m2 die Hevrer bringt, sind jedoch infolge einer 
7 ungliicklichen Auswahl der Parameter ohne 
Wert. AuBerdem wurde in dieser Arbeit kein 


Hoy Z b) Unterschied zwischen evolutionaren und nicht- 
thy evolutionéren StoBwellen gemacht. 

2 Tn der spater erschienenen Arbeit von Ltst [46] 

7 A ae wird eine Untersuchung der Parameterabhan- 

Lokee gigkeit wesentlicher physikalischer GroBen 

ons sie der StoBwelle (Dichte, Druck, Magnetfeld, 

igs *, 3a Warmeenergie und Winkel zwischen der Rich- 

aire 


WWE tung des Magnetfeldes und der Normalen) 
4 8 120 16h, ss : DN rie tae “ 
durchgeftihrt. Die Abhangigkeit dieser GroBen 
Bild 3. Dichte (a) und magnetischer Feld- : A . a Il F 
stirke (b) in der StoBwelle in Ab- Wurde in vielen Diagrammen dargestellt, die 
hiingigkeit von der Mach-Zahl. Die mit einer elektronischen Rechenmaschine er- 
nichtevolutioniren Teile sind punk- : a x ; 
fiery peneichiel, (uncut eiaiates halten wurden. Leider berticksichtigt der Autor 
die schnelle, C_ die langsame Sto8- die Evolutionsbedingungen nicht; in Wirk- 
welle; C4 markiert die (nicht wirk- Jichkeit kénnen nur die Abschnitte der von 
lich existierende) nichtevolutionire x rhe Aas j 
Stopwolle. Lust aufgestellten Kurve existieren, die evo- 
lutionaren StoBwellen entsprechen. Als Bei- 
spiel sind auf Bild3a und b die von List gefundenen Abhangigkeiten der 
GroBen 02/0, und H,,/H,, von der Mach-Zahl 


2 


v 
Mee ies 
=A U2 
dargestellt; dabei ist 
3p,/2 dake 


(34/2) + (7/82) 


Kine vollstandigere Untersuchung der GréBen in der magnetohydrodynamischen 
StoBwelle wurde von BazEr und Ericson [47] durchgefiihrt. Jedoch auch in 
dieser Arbeit finden die Evolutionsbedingungen keine Beriicksichtigung. Im 
folgenden werden wir die Resultate von BazEr und Ericson durch Untersuchun- 
gen evolutionarer Teile der StoBadiabate erganzen [55]. 

Ks gibt schnelle und langsame StoBwellen. Schnelle Wellen sind immer evolu- 
tionar. Langsame StoBwellen von kleiner Intensitaét sind auch evolutionar. In 
dem MaBe, wie sich der Dichtesprung Ao = 0, — 0, der langsamen StoBwelle 
vergroBert, wird das transversale Magnetfeld H ay hinter der Welle geschwacht, 
und in der Grenze des Evolutionsgebietes geht es gegen Null. Bei weiterem Fort- 
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schreiten lings der StoBadiabate folgt der nichtevolutionare Teil der langsamen 
StoBwellen, der stetig zum Alfvénschen Sprung fiihrt (siehe Bild 2). Druck- und 
Entropiesprung in schnellen und langsamen StoSwellen sind im Evolutionsgebiet 
monoton wachsende Funktionen des Dichtesprunges. In schnellen StoBwellen 
ist der maximale Dichtesprung gleich (A @)max = 20,/(y — 1), wobei » der 
Poissonsche Adiabatenexponent ist; hierbei geht der Druck- und der Entropie- 
sprung gegen unendlich. Es gibt zwei Typen der Abhangigkeit des Magnetfeld- 
sprunges 4H, = H,, — H,, vom Dichtesprung in schnellen Wellen. In Wellen 
des ersten Typs, die existieren, wenn : 


sin2 6, > (y o 1) (1 , r,) 
if 
ist, (ry = c{/ Uj, = 42yp,/Hz, 0,: Winkel zwischen der Richtung des Magnet- 


feldes H, und der Normalen der Sprungflache) wachst der Magnetfeldsprung beim 
Anwachsen des Dichtesprunges monoton von Null bis zum maximalen Wert 


2H. 
(4H) max — a 
an. 
In Wellen des zweiten Typs, die existieren, wenn 
a4 | HG ee 
sin2 6, zx (y ) ( r,) 


a 


ist, gibt es eine nichtmonotone Abhangigkeit des Magnetfeldsprunges vom 
Dichtesprung: bei Anwachsen des Dichtesprunges vergréBert sich anfangs der 
Magnetfeldsprung von Null bis zu einer gewissen MaximalgréBe, und danach 
sinkt er wieder bis auf den Wert 


Auf dem evolutionéren Abschnitt der langsamen StoBwelle wachst der Feld- 
sprung immer mit dem Dichtesprung?). 

Im Grenzfalle 6, > 0,7, > 1 verhalt sich die schnelle Welle so, als sei das Magnet- 
~ feld nicht vorhanden, und die langsame Welle hat eine unendlich kleine Amplitude. 
Im Falle von 6, +0, 7; <1 gehdrt die schnelle Welle zum zweiten Typ. Bei 
kleiner Intensitét der schnellen StoBwelle (92/01 < (y + 1 —2r,)/(y — 1)) ist das 
transversale Magnetfeld H,, hinter der Welle noch von Null verschieden. Wenn die 
Intensitat der StoBwelle o,/0, den Wert (y + 1 — 2r,)/(y — 1) tibersteigt, wird 
das transversale Magnetfeld H,, hinter der Welle gleich Null, und die StoBwelle 
verhalt sich so wie beim Nichtvorhandensein des Magnetfeldes. 

Die langsame Welle verhalt sich im ganzen evolutionaren Bereich fur 0 i> 9, 
r, <1 so, als sei kein Magnetfeld vorhanden. Im Grenzfalle 6, > 2/2 ist die 
schnelle StoBwelle einer Welle des ersten Typs, und die langsame Stofwelle ver- 
wandelt sich in einen tangentiellen Sprung. 


1) Mit anderen Worten: Langsame StoBwellen gehoren immer zum ersten Typ. Die Behauptung 
von Bazer und Ericson [47], in langsamen Stofwellen sei eine nichtmonotone Abhangigkeit 
des Magnetfeldsprunges vom Dichtesprung méglich, wurde ohne Beriicksichtigung der Evolu- 
tionsbedingungen aufgestellt. 
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Das Vorhandensein eines Magnetfeldes vergréBert den Drucksprung bei fest vor- 

gebenem Dichtesprung [47]?). ' 

Tn schnellen StoBwellen groBer Intensitaét (p, > p, + H 2/8 7) ist das Vorhanden- 

sein eines Magnetfeldes unwesentlich [53, 46]. Insbesondere ist die gréBte Kom- 
pression 0,/0,, die in der StoBwelle erhalten 

a) werden kann [57, 47], gleich (y + 1)/(y —1). 

Pe e Eine langsame StoBwelle kann keine beliebig 

groBe Intensitat haben. Deshalb gibt es bei 


Pp. > Pp, + Hi/8a 


nur eine StoBwelle [36, 46] (naémlich eine 
schnelle). 
In den Veranderlichen (1/0, 72) wurden evo- 
lutiondre Bereiche der StoBadiabate nur fiir 
die Grenzfalle H,,<H, (Bild4a) und H,, >H, 
(Bild 4b) bestimmt [55]. Wie aus Bild 4a er- 
sichtlich ist, ist im Falle von H,,< H, der 
@ b) Abschnitt der StoBadiabate, der einer schnellen 
StoBwelle entspricht, hdher gelegen, als der 
Bereich, der einer langsamen StoBwelle ent- 
spricht. Das bedeutet, daB bei H,,< H, die 
schnelle StoBwelle thermodynamisch giinsti- 
ger ist: auf ihr wachst die Entropie sehr 
stark. Im Falle von H,, >> H, erweist sich 
die langsame StoBwelle als thermodynamisch 
ile ihe all giinstiger (Bild 4b). 
40, f, Pp Den Unterschied in der relativen Lage der 
Bild4. Die evolutionaren Teile der Stofadiabate 'Teile der Adiabate, die einer schnellen und 
Ne ae hee einer langsamen StoBwelle entsprechen, kann 
kennzeichnet die schnelle, das Zeichen — man auf Bild 4a und 4b auf folgende Weise 
oe ame Sobel. erklaren. In der StoBwelle geht ein Uber- 
gang der kinetischen, Energie 0,v?,/2 in ma- 
gnetische H3/8x und Warmeenergie 3p,/2 vor sich. Der schnellen StoBwelle 
entspricht eine gréBere Ubertragung an kinetischer Energie. Im Falle von 
H,,< H, wird die magnetische Energie wenig geiindert; deshalb entspricht der 
schnellen StoBwelle eine gréBere Aufheizung. Im Falle von H,, > H, wird die ma- 
gnetische Energie merklich beeinflu8t. Da das Magnetfeld in einer schnellen 
StoBwelle verstarkt wird, in einer langsamen aber geschwacht wird [38], ist es klar, 
dafs in einer schnellen StoBwelle bei H,, > H, eine Aufheizung vor sich geht. 


IN 
le 
> 
Shs 


§ 7 Das Kolbenproblem 


Bei Betrachtung des Kolbenproblems tauchen in der Magnetohydrodynamik zwei 
prinzipielle Schwierigkeiten auf, die es in der gewohnlichen Hydrodynamik nicht 
gibt. Die erste Schwierigkeit besteht darin, daB die Grenzbedingungen an den 
Sprungflachen und auf dem Kolben die entstehenden StoBwellen und stetigen 
Wellen nicht eindeutig bestimmen. Im Falle eines idealen Gases erhalt man eine 


') Im Falle H, = 0 wurde diese Tatsache von G. 8S. Goxtcyn [56] festgestellt. 
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groBe Zahl von Lésungen fiir eine vorgegebene Bewegung des Kolbens. Es ent- 
spricht jedoch nur eine dieser Lésungen den Evolutionsbedingungen. Deshalb 
beseitigt die Elimination der nichtevolutiondren Spriinge die oben erwahnte 
Schwierigkeit und lést das Kolbenproblem mathematisch korrekt.") 

Kine zweite Schwierigkeit besteht darin, daB 
stetige Strémungen, die bei der Bewegung des 
Kolbens entstehen, nicht immer einfache Wel- 
len sind. Selbstverstandlich findet man beim 
Fehlen von Spriingen an eine konstante Stré- 
mung angrenzend zunachst einfache Wellen. 
Jedoch braucht dieses Gebiet nicht bis zum 
Kolben zu reichen, denn zwischen dem Kolben 
und der einfachen Welle kann ein zweites 
Gebiet liegen, in dem die Welle nicht einfach 
ist (Bild 5). Um diese Schwierigkeit. zu um- 
gehen, werden wir nur Bewegungen des Kol- 
bens mit konstanter Geschwindigkeit unter- 
- suchen. Da hierbei im Problem ein Parameter 
fehlt, der die Dimension einer Lange hat, wird 
die Bewegung des Mediums automodellartig 
sein, d. h., alle entstehenden Wellen werden 
entweder Sprungwellen oder einfache Wellen 
sein. 

Die Bewegung des Mediums unter dem EinfluB 
der Kolbenbewegung wird charakterisiert durch 


Ch 


b) x 


Bild 5. Wellen, die bei der Bewegung eines 


die Reihenfolge der nacheinanderlaufenden 
~StoB- und Automodellwellen. Wie wir gesehen 
haben, gibt es drei Arten von evolutionaren 
Sprungwellen: schnelle und langsame StoB- 
wellen und Alfvénsche Spriinge. AuBerdem 
existieren zwei Arten von stetigen Lésungen: 
schnelleund langsame magnetoakustische Auto- 
modellwellen. Wie weiter oben bemerkt wurde, 
sind StoBwellen Kompressionswellen und Auto- 
modellwellen Dilatationswellen. Die Geschwin- 
digkeiten dieser Wellen sind so, daB vorn eine 
schnelle Welle laufen kann (StoB- oder Auto- 


Kolbens auftreten: 

a) in der gewodhnlichen Hydrodyna- 
mik, b) in der Magnetohydrodyna- 
mik. K: Bahnkuryve des Kolbens, 
Ch: Charakteristik in der gewohn- 
lichen Hydrodynamik da/dt=v,+¢; 
Ch,.: ,,5chnelle* Charakteristik in 
der Magnetohydrodynamik da/dt = 
Vz + U,; Ch_: ,,Langsame“ Charak- 
teristik in der Magnetohydrodyna- 
mik da/dt = vz; + U_; 0: Ungestértes 
Gebiet; 1: Hinfache Welle; 2: Welle 
allgemeineren Typs (nicht mehr ein- 
fach). 


modellwelle), ihr folgt ein Alfvénscher Sprung und schlieBlich eine langsame Welle 
(StoB- oder Automodellwelle). Wenn man beriicksichtigt, daB einige der aufgezahl- 
ten Wellen fehlen kénnen, so erhalt man eine groBe Zahl qualitativ verschiedener 
Bewegungsbilder des Mediums, die bei verschiedenen Werten der Sprunggeschwin- 
digkeit realisiert werden. 
Auf der Kolbenoberflache gelten offensichtlich die Grenzbedingungen 

Vg = Ug: (7.1) 


Um die restlichen zwei Bedingungen zu erhalten, mu8 man zu einem Bezugssystem 
iibergehen, das sich mit dem Kolben mitbewegt, wobei auf der Oberflaiche des 


1) Diese Schwierigkeit kann man auch ohne Benutzung der Evolutionsbedingungen umgehen, 
e : ; = 
wenn man sich auf Kolbenbewegungen beschrankt, bei denen keine langsamen StoBwellen 


entstehen [58]. 
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aaa : ; = — Iten (der Strich 
toefiihrten Kolbens die Grenzbedingungen FH, 0, Ey = 0. ge | 
ae daB die entsprechende Gré8e in dem mit dem Kolben mitbewegten 
Koordinatensystem gemessen ist). Infolge der unendlichen Leitfaihigkeit des 
Mediums ist Z’ = — 1/c[v’, H] (c: Lichtgeschwindigkeit), und daraus folgt 


0, =u, (H, #9). 


Shek 4 : niakarenes 
Aut diese Weise ist auf metallischen Oberflichen die Relativgeschwindigkeit 

leitenden Flissigkeit gleich Null (eine Ausnahme bildet der Fall, daB das Magnetfeld 
keine Normalkomponente hat). Es zeigt 


(7.2) 


Vy = Uy, 


man _,,elektrodynamischeZihigkeit [59] 
nennt. 

Bei geniigend groBer Amplitude einer 
DilatationsstoBwelle geht die Dichte 
des Mediums hinter der Welle gegen 
Null, es tritt Kavitation ein. Hierbei 
miissen an der Grenze des Vakuums 
Magnetfeld und transversale Kompo- 
nente des elektrischen Feldes stetig sein, 
d.h., sie miissen die Grenzbedingungen 


e = 0, 
HH, (uy, — vy) — 2, a Se 
iif (uz je Vz) 7. H, (Uz Vz) =0 (7.3) 
erftillen. 


Betrachten wir nun den interessanten 
Fall, daB das Magnetfeld, die Kolben- 
geschwindigkeit und die Normale der 
Kolbenoberflache in einer Ebene liegen 


Bild 6. Wellen, die bei der Bewegung eines Kolbens 
entstehen. 
Auf der Abszisse ist die Lingskomponente der 
Kolbengeschwindigkeit u,, auf der Ordinate 
die Querkomponente w, aufgetragen. Die Buch- 
staben Ut, U-, R*+, R-, A kennzeichnen das 
Vorhandensein von schnellen und langsamen 
StoBwellen, schnellen und langsamen Dilata- 
tionswellen (Automodellwellen) bzw. die Exi- 
stenz eines Alfvénschen Sprunges. GV ist der 
Punkt gr6Bter, in der Automodellwelle erreich- 
ter Dilatation. W kennzeichnet das Auftreten 
von Kavitation. 


(dies sei die xy-Ebene). Die GroBen v, 
und H, werden hierbei nicht nur im 
nichtbewegten Medium, sondern auch 
in allen entstehenden Wellen gleich 
Null sein. Deshalb kann der Alfvénsche 
Sprung das Magnetfeld nur um 180° 
drehen. Machen wir noch die ungestérte 
Geschwindigkeit des Mediums v, zu 
Null und setzen wir auBerdem fest, 


sich eine interessante Erscheinung, die — 


<< 


da die Komponenten des nichtgestér- 
ten Magnetfeldes Hy, und Ho, positiv 
sind, so hangen die Wellen, die bei der Bewegung des Kolbens unter diesen 
Bedingungen entstehen, nur noch von seiner Geschwindigkeit (w,, Uy) ab. Diese 
Abhangigkeit [60] ist in Bild 6 dargestellt. Bei einer gentgend groBen Am- 
plitude der langsamen Verdiinnungswelle geht die Dichte des Mediums hinter 
der Welle gegen Null, man erhilt Kavitation. Im Vergleich zur gewohnlichen 
Hydrodynamik, in der die Kavitation eintritt beim Weegziehen des Kolbens mit 
einer Geschwindigkeit, die 2c /(y — 1) wtihbersteigt (cy: Schallgeschwindigkeit im 
nichtgestorten Medium), entsteht in der Magnetohydrodynamik die Kavitation 
auch bei geringeren Kolbengeschwindigkeiten, wenn nur die Bewegungsgeschwin- 
digkeit des Kolbens in Querrichtung gentigend gro8 ist. Wenn sich der Kolben nur 
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in Querrichtung bewegt, so tritt die Kavitation ein, wenn die Kolbengeschwindig- 
keit die Schallgeschwindigkeit im ungestérten Medium [29] bei y = 5/3 um den 
Faktor 3,67 lbersteigt. Die Kavitation entsteht auch, wenn der Kolben in das 
Medium hineinbewegt und gleichzeitig in Querrichtung bewegt wird. Beim Ein- 
schieBen des Kolbens mit Uberschallgeschwindigkeit [60] tritt die Kavitation ein, 
wenn der Winkel zwischen dem Geschwindigkeitsvektor des Kolbens und der Nor- 
male seiner Oberflache groBer wird als 70° (fiir y = 5/3). In diesem Falle ist die 
Differenz zwischen der Geschwindigkeit, mit der sich die Grenze Medium —Vakuum 
verschiebt, und der Kolbengeschwindigkeit sehr klein. 
Im Gegensatz zur langsamen Verdiin- 
nungswelle ist in einer schnellen Ver- H 
diinnungswelle Kavitation unméglich. f 
Wenn im nichtgestérten Medium die uy 
Alfvénsche Geschwindigkeit wesentlich a Rave 
kleiner ist als die Schallgeschwindigkeit, | y U 
so wird bei maximaler Amplitude der a) 
schnellen Verdiinnungswelle die Alfvén- 
sche Geschwindigkeit hinter ihr von ahn- 
licher GréBe sein wie die Schallgeschwin- 
digkeit. Bei EKinhaltung der Bedingung 
H?/8a<p, wird die Dichte des Me- “y 
diums im Punkte der gré8ten Verdiin- 
nung (Punkt GV auf Bild 6) auBerst 
klein sein. 
Bei Uberschallgeschwindigkeit des Kol- C) d) 
bens tritt Erzeugung eines Magnetfeldes Bild 7. Ausbildung von StoBwellen (U), Alfvénschen 
ein, d.h. ein Anwachsen des Magnet- PSE NT cor etd bre aad 
t ax $ pannungen bei Quer 
feldes von unendlich kleiner zu end- bewegung des Kolbens, 
licher GroBe ; hierbei kommt der magne- 
tische Druck in die GréBenordnung des hydrostatischen oder tibertrifft ihn. Beim 
Hineinschieben und Gleiten des Kolbens mit Uberschallgeschwindigkeit ist das 
erzeugte Magnetfeld der Normalkomponente der Kolbengeschwindigkeit direkt 
proportional. 
Die topologische Struktur des Bildes 6 kann man auch ohne Rechnung auf Grund 
qualitativer Uberlegungen erhalten [61]. Dazu legen wir fest, durch welche Gebiete 
die Ordinatenachse u,—0 verlauft. Da die magnetischen Kraftlinien an das 
Medium und an den Kolben ,,angefroren“ sind, werden sie deformiert. Fir u, < 0 
zeigt dies Bild 7a. Die Kriimmung der magnetischen Kraftlinien fihrt zur Er- 
scheinung quasielastischer Spannungen, die die Kraftlinien zu glatten versuchen 
(siehe den Pfeil auf Bild 7a). 
Weil in der Nahe des Kolbens und im Unendlichen vy, = 0 ist, entsteht vor dem 
Pfeil eine Kompressionswelle (StoBwelle) und dahinter eine Verdiinnungswelle 
(Automodellwelle). Eine Alfvénsche Welle entsteht im vorliegenden Falle nicht, 
weil H, in der Nahe des Kolbens und im Unendlichen gleiches Vorzeichen hat. 
(Erinnern wir uns, da8 in Sto8- und Automodellwellen das Vorzeichen von H, 
nicht gedndert wird [38]; beim 180°-Alfvén-Sprung aber andert es sich.) Auf diese 
Weise liuft bei w, = 0, uy < 0 vorn eine StoBwelle, hinter ihr eine Automodellwelle 
(Kombination Ut+R-). 
Analoge Uberlegungen zeigen, daB in dem auf Bild 7b dargestellten Falle die 
Kombination R+U- auftritt. Bei VergréBerung der Geschwindigkeit uw, andert 


. 
U \ 
A R 
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sich das Vorzeichen von H, in der Nahe des Kolbens. Das zieht das Auftreten einer 
Alfvénschen Welle (R+AU-, Bild 7c) nach sich. Bei weiterer VergroBerung Von Uy 
wird die GréBe H, in der Nahe des Kolbens gréSer als der Wert H, im Unend- 
lichen; dann ist die Resultierende der Spannungen vom Kolben weggerichtet 
(Bild 7d), das entspricht der Kombination U+tAR-. Bei gentigend groBem Wert 
von u, wird die Amplitude der Dilatationswelle so groB, dals die Dichte des 
Mediums hinter der Welle-gleich Null wird; es tritt Kavitation ein. — ; 
Ein Bild der Bewegung des Mediums bei uv, 4 0, uy ~ 0, u, = 0 kann man erhalten, 
wenn man vom im vorstehenden betrachteten Falle u, = 0 ausgeht. Bei einer — 
VergroBerung von wu, wird die Amplitude der Dilatationswelle kleiner und die 
Amplitude der Kompressionswelle groBer1). Bei einem gewissen Wert von uw, ver- 
wandelt sich die Dilatationswelle in eine Kompressionswelle. Analog fihrt eine 
Verkleinerung von u, zur Wandlung der Kompressionswelle in eine Dilatations- 
welle. Hine weitere Verkleinerung der GroBe wu, fiihrt zum Auftreten der Kavi- 
tation. Bei Bewegung des Kolbens tritt eine Bremskraft F auf, die zwei Kompo- 
nenten hat; die Langskomponente F’, (Frontwiderstand) und die Querkomponente 
F,, (Schubkraft). Bei einer schnellen StoBwelle mit groBer Amplitude, die im 


Bereich 
i > Cy; | uy | << Uz, Hj/8a K< Po 


von einer langsamen Dilatationswelle begleitet wird, wird der Frontwiderstand 
durch die Formel 


_ (y +1) eur 
7 2G 1) 
bestimmt [60]. 
Wenn man | w,| bei fixiertem wu, vergréBert, so wird die GréBe des Frontwiderstan- 
des kleiner und nimmt auf der Kavitationslinie den Wert 


Be ee oS. 


an, 


Bei weiterer VergréRerung von | w,| verandert sich die GréBe des Frontwiderstandes 
nicht. 


Die Schubkraft wird bestimmt durch die Formeln 


i he (y+ sun aes 
ve 


mit | uy| << uw, 


bzw. im Falle der Kavitation 


et Vy+1 009 Vor Sign uy. 


Der Fall, in dem die Geschwindigkeit des Kolbens nicht in der Ebene des Magnet- 
feldes (xy) liegt, d.h., wenn wu, 4 0 ist, unterscheidet sich vom vorhergehenden 


dadurch, daB die Richtungsinderung des Magnetfeldes beim Alfvénschen Sprung 
von 180° verschieden ist. 


*) Eine Ausnahme von dieser Regel findet man im Bereich U+U- bei Uberschallgeschwindig- 
keit des Kolbens. Bei einer VergréBerung von u x geht eine Neuverteilung der Amplituden von 
schneller und langsamer Sto8welle vor sich; das Anwachsen der Amplitude der schnellen 
StoBwelle ist dabei begleitet von einer Verringerung der Amplitude der langsamen StoBwelle. 
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§ 8 Zerfall eines Sprunges 


Es erhebt sich die Frage: Was geschieht mit einem Sprung in den Anfangsbedin- 
gungen, wenn auf ihm die notwendigen Grenzbedingungen (Kontinuitat der 
Massenstr6mung, des Impulses usw.) nicht erfiillbar sind? Solche Anfangsbedin- 
gungen erhalt man z. B. beim Zusammenprall von Gasmassen oder bei plotzlicher 
Zerstérung einer Scheidewand zwischen zwei Gasen, die unter verschiedenem 
Druck stehen. Dieses Problem wurde von RrmMaNnn aufgeworfen und fiir die 
gewohnliche Hydrodynamik von N. E. Korcurnn |62] gelést. Wie die Unter- 
suchung zeigt, spaltet der Sprung in drei Wellen auf [34]: eine Welle, die sich nach 
rechts bewegt (Sto8- oder Automodellwelle), eine Welle, die sich nach links bewegt 
(StoB- oder Automodellwelle) und einen tangentialen Sprung. 


In der Magnetohydrodynamik verhalt sich die Lésung des oben genannten Pro- 
blems beim Fehlen eines longitudinalen Magnetfeldes H,, qualitativ ebenso wie 
beim Fehlen eines Magnetfeldes tiberhaupt [63]. Der Zusammenprall zweier Gas- 
massen mit H, = 0 und nachfolgendem Zerfall des Sprunges wurde von T. F. Vou- 
KOV [64] im Zusammenhang mit dem Problem der Plasmaaufheizung bis zu thermo- 
nuklearen Temperaturen durch StoBwellen behandelt. 


Der Zerfall des Sprunges im ruhenden Plasma mit H, = 0, |Ap*| < p* (p* = p 
+ H?/8x) wurde von Kato [65] untersucht. Die erhaltene Lésung wurde zur 
qualitativen Betrachtung der Pulsation eines Plasmazylinders im longitudinalen 
Magnetfeld verwendet. 


Wenn das longitudinale Magnetfeld H, von Null verschieden ist, so ist die Spal- 
tung des Ausgangssprunges in der Magnetohydrodynamik von der Aufspaltung 
in der gewohnlichen Hydrodynamik qualitativ verschieden. Der Anfangssprung 
spaltet jetzt in sieben Wellen auf: drei Wellen, die sich nach rechts bewegen, drei 
Wellen, die sich nach links bewegen und zwischen ihnen ein Kontaktsprung, der 
relativ zum Medium ruht. 

Jede dieser Wellen wird durch einen Parameter charakterisiert (durch eine 
Amplitude, d. h. durch den Sprung einer der magnetohydrodynamischen GréBen). 
Andererseits ist die Summe aller Spriinge einer jeden magnetohydrodynamischen 
GroBe in den sieben entstehenden Wellen gleich dem Anfangssprung. Weil die 
Zahl der magnetohydrodynamischen GréBen gleich sieben ist (0, p, Vz, Vy, Vz, Hy, 
-H,), so erhalten wir ein System von sieben Gleichungen mit sieben Unbekannten, 
bei dessen Lésung wir die Amplituden aller Wellen finden, die beim Zerfall des 
Anfangssprunges entstehen. 

Die allgemeine Lésung des Zerfallsproblems wurde wegen groBer mathemati- 
scher Schwiergikeiten nicht gefunden. Fiir den Fall, da& der Anfangssprung sehr 
klein ist!), wurde das Problem gelést [66]. In diesem Falle sind alle sich bildenden 
Spriinge ebenfalls klein. Da das Sprungverhaltnis der magnetohydrodynamischen 
GréBen in einer Automodell- und einer StoBwelle kleiner Intensitat in erster Nahe- 
rung gleich ist, besteht der Unterschied zwischen einer StoBwelle und einer Auto- 
modellwelle nur darin, daB in der StoBwelle die Dichte wachst, wahrend sie in 
der Automodellwelle abnimmt. Die Dichtespriinge in schnellen und langsamen 


1) Der allgemeinere Fall, daB das System der magnetohydrodynamischen Grundgleichung 
durch ein hyperbolisches System aus n Gleichungen mit 7 Unbekannten ersetzt wird, wurde 


von Lax [19] untersucht. 


aie 
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magnetoakustischen Wellen werden durch die Formel 
ec? UU; | do — ge As : f 
dines Ay , @Uz(HiAu , Uidvs 
pe (orto Heel eee 
(8.1) 


) 
As ar oR (Siete Sarge 82 U. 


gegeben, dabei ist R= y (U2 + c?)? — 4c? Uz, Ao: Dichtesprung, 4s: Entropie- 
sprung, Av: Geschwindigkeitssprung und A H;: Sprung des transversalen Magnet- 
feldes im Anfangssprung; die Zeichen + haben denselben Sinn wie in (2.2). Mit 
der angefiihrten Formel kann man die Vorzeichen von ae o bestimmen, und 
damit erklaren, in welche Wellen der Anfangssprung zerfallt. Der StoBwelle ent- 


spricht ein Ae o > 0 und der Automodellwelle ein ae oa 0. . 

Der Sprung in den Anfangsbedingungen hat Ahnlichkeit mit einer nichtevolu- 
tiondren StoBwelle. Obgleich im letzteren Falle alle Grenzbedingungen erfillt 
werden, gentigt eine unendlich kleine St6rung, um die Grenzbedingungen zu 
verletzen. Hine Verletzung der Grenzbedingungen fiihrt zu einer Spaltung der 
nichtevolutionaren StoBwelle in einige auslaufende Wellen (Sprungwellen oder 
stetige Wellen). Hine solche Spaltung wurde in [45] untersucht fiir den Fall, daB 
das Magnetfeld an beiden Seiten der StoBwelle einen kleinen Winkel mit der Nor- 
male der Sprungfliche bildet und die Ausbreitungsgeschwindigkeit der StoB8welle 
Vy, nahezu der Alfvénschen Geschwindigkeit U,, entspricht, die gréBer ist als die 
Schallgeschwindigkeit c,. (Der Index 1 bezieht sich auf das Gebiet vor der Welle.) 
Die Evolutionsbedingung einer langsamen StoBwelle-ist in diesem Falle: 71, < Uj,. 
Bei einer Storung dieser Bedingung wird die StoBwelle nichtevolutionar. Eine 
solche nichtevolutionare StoBwelle kann man erhalten, wenn man eine gewohn- 
liche stabile hydrodynamische StoBwelle, bei der v,, > c, und vy, < Cy ist, in ein 
Magnetfeld bringt. Wenn das transversale Magnetfeld H, gleich Null ist, so kann 
die betrachtete StoBwelle in zwei getrennte StoBwellen gespalten werden [42], 
wobei die magnetohydrodynamischen Groen in dem zwischen diesen Wellen 
eingeschlossenen Gebiet gegeben sind durch: 


5 G1 rape ls 
Ui, ‘ Vie 
a ee y2 2 oe 2 
b =p, + Se = Uae) = + Me = Vis) 
12 
8.2 
5 4 [2lote — Ute) (4U2, — ob, — 3e8) Se 
je as 820, (vi, = Uir) (4 Cee — mn a 3c) Ia 
ier SUE, ata 


Kine solche Spaltung ist nur méglich [45], wenn die Evolutionsbedingung v,, < U,, 


nicht erfullt ist. Das folgt daraus, daB bei v,, < U yr die Ausdriicke fiir v, und H 
imaginar werden. d : 


Beriicksichtigt man ein schwaches transversales Magnetfeld, so wird die nicht- 
evolutionére AusgangsstoBwelle in vier Wellen gespalten (wenn man Wellen 


oe 
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nicht beriicksichtigt, deren Amplituden die Gré8enordnung a«H,, haben, 


6= V (Ure — Uy,y)/Uy, <1). Man erhalt eine nach links (die w-Achse zeigt nach 
rechts) laufende schnelle StoBwelle mit einer Amplitude der GréBenordnung « 
und einer Geschwindigkeit der Gré8enordnung Uy, = H,,/ 4201; eine nach links 
laufende langsame StoBwelle mit einer Amplitude, die sich wenig von der Ampli- 
tude der Ausgangswelle unterscheidet und mit einer Geschwindigkeit der GréBen- 
ordnung « H,, lauft; auBerdem einen nach rechts laufenden Alfvénschen Sprung, 
der das Magnetfeld um 180° dreht, und schlieBlich eine nach rechts laufende 
schnelle StoBwelle mit einer Amplitude der GréBenordnung H 
Mit dem Zerfall des Sprunges in den Anfangs- 

bedingungen ist die Frage der Uberginge 4,49 UtA U7 K UA Ut 


zwischen magnetohydrodynamischen Spriin- 

gen eng verbunden?). Die Beriicksichtigung 

der Evolutionsbedingungen verindert ein 

wenig das Bild der méglichen Ubergange. x Es 
H.=0 U* T U* 


Vor allem erweist sich die Behauptung, daB 
der Alfvénsche Sprung stetig in eine StoBwelle 


ly: 


Bild 8. Schema der Ubergiinge zwischen ma- 


tibergehen kann, als unwahr. Denn der Alfvén- gnetohydrodynamischen Spriingen. 
oe . . i F a = 
sche Sprung kénnte mit einer StoBwelle nur Die Butane eee 
kennzeichnen schnelle und langsame 
dann zusammenfallen, wenn das Magnetfeld StoBwellen, Alfvénsche, Kontakt- und 
bis zum Sprung und nach ihm in einer Ebene Tangenualepringe- 


liegt. Im Alfvénschen Sprung wird dann das 

Magnetfeld um 180° gedreht. Aber dabei verandert die transversale Kompo- 
nente des Magnetfeldes das Vorzeichen, wahrend sich auf einer StoBwelle 
das Vorzeichen der transversalen Komponente des Magnetfeldes nicht andert. 
Deshalb ist ein Ubergang zwischen einem Alfvénschen Sprung und einer evolu- 
tionaren StoBwelle nicht méglich. Es geht aber eine (nicht wirklich existierende) 
nichtevolutionare langsame StoBwelle formal stetig in einen Alfvénschen Sprung 
tiber [47] (siehe Bild 2). 

Zwischen einer schnellen und einer langsamen StoBwelle sind stetige Ubergiinge 
ebenfalls unméglich. Das folgt daraus, daB die Bereiche, in denen schnelle und 
langsame StoBwellen existieren, nicht zuasammenhangen (siche Bild 2). 

Ebenso kann eine schnelle StoBwelle nicht stetig in einen tangentiellen Sprung 
iibergehen, weil das den Bedingungen (4.3) widersprechen wiirde. Deshalb sind 
-Ubergiinge nur zwischen tangentiellen und Kontaktspriingen, zwischen tangen- 
tiellen und Alfvénschen Spriingen und zwischen einem tangentiellen Sprung und 
einer langsamen StoBwelle méglich. Die Bedeutung méglicher Ubergiinge zwischen 
magnetohydrodynamischen Spriingen wird klarer, wenn man das Problem des 
Zerfalls willkiirlicher Spriinge in den Anfangsbedingungen untersucht. Wenn das 
normale Magnetfeld H,, von Null verschieden ist, so zerfillt der Sprung in sieben 
Wellen, jede von ihnen wird durch einen Parameter charakterisiert. 

Wenn das normale Magnetfeld /, gleich Null ist, so zerfallt der Anfangssprung in 
drei Wellen: eine schnelle StoBwelle, die sich nach rechts bewegt; eine schnelle 
StoBwelle, die sich nach links bewegt, und einen tangentiellen Sprung zwischen 
ihnen. Jede der StoBwellen wird durch einen Parameter charakterisiert, der 
tangentielle Sprung durch fiinf Parameter [30]. Die Gesamtzahl der Parameter 
ist gleich sieben, d.h. gleich der Zahl der Spriinge magnetohydrodynamischer 


GréoBen auf dem Anfangssprung. 


1) Diese Frage wurde erstmals in der Arbeit [32] von S. I. SyrovatsxK1s behandelt. 
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Aus diesem Grunde’ stellt der tangentielle Sprung die Vereinigung von finf 
Spriingen dar (zwei langsame StoBwellen, zwei Alfvénsche Spriinge, ein Kontakt- 
sprung). Ein Schema der Uberginge zwischen magnetohydrodynamischen 
Spriingen ist auf Bild 8 dargestellt. 
Ubersetzt von L. u. M. RorHHarptT 
Deutsche Red. durch L. RoTHHARDT 
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